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Abstract 

In this work, we perform the one-loop calculation of the scalar and pseudoscalar form factors in the frame- 
work of U(3) chiral perturbation theory with explicit tree level exchanges of resonances. The meson-meson 
scattering calculation from Ref. [1] is extended as well. The spectral functions of the nonet scalar-scalar (SS) 
and pseudoscalar-pseudoscalar (PP) correlators are constructed by using the corresponding form factors. Af- 
ter fitting the unknown parameters to the scattering data, we discuss the resonance content of the resulting 
scattering amplitudes. We also study spectral-function sum rules in the SS-SS, PP-PP and SS-PP sectors 
as well as semilocal duality from scattering. The former relate the scalar and pseudoscalar spectra between 
themselves while the latter mainly connects the scalar spectrum with the vector one. Finally we investigate 
these items as a function of Nc for Nc > 3. All these results pose strong constraints on the scalar dynamics 
and spectroscopy that are discussed. They are successfully fulfilled by our meson-meson scattering amplitudes 
and spectral functions. 
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1 Introduction 



Spectral functions of current-current correlators are interesting objects in hadron physics. They are sensitive to 
both low and high energy dynamics of QCD, which can be evaluated in two reliable approaches: chiral perturbation 
theory (xPT) [2, 3, 4] in the low-energy region and the operator product expansion (OPE) [5, 6] for the higher 
one. The celebrated Weinberg sum rules are in fact derived by studying the differences of the vector-vector 
and axial- vector-axial- vector spectral functions [7], where not only the low and asymptotic energy regions but 
also the resonance region are considered. After then, great progresses, both on the experimental and theoretical 
sides [8], have been made along this research line. The vector and axial- vector spectral functions are experimentally 
measurable quantities, mainly through the r decays, due to the noticeable V — A nature of the Standard Model in 
the energy region well below the W boson mass. 

Concerning the scalar and pseudoscalar spectral functions, there are no direct experimental data. However, 
great efforts have been made on the theoretical sides both in perturbative QCD calculations [9, 10, 11], which 
arc important to reduce the QCD background in the search of the Higgs particle, and in the nonperturbative 
QCD region [3, 12, 13, 14, 15, 16, 17, 18, 19, 20]. Through the matching between the low and high energy 
behaviors of the spectral functions, valuable information can be obtained: light quark masses and Cabibbo angle 
\V US \ [9, 18, 19], determination of the low energy constants (LECs) in xPT [12, 13, 16, 17, 21] and constraints on 
resonance parameters [14, 16, 17, 20, 22, 23]. 

On the other hand, in the past decade scalar dynamics in the nonperturbative QCD region has been greatly put 
forward thanks to the combination of x?T and nonperturbative approaches from the S'-matrix theory [24, 25, 26, 
27, 28, 29, 30, 31, 32, 33, 34, 35]. One of the main results is the reappearance of the broad /o(500) (traditionally 
called a), the lightest resonance in QCD spectrum, after its long absence in the PDG list [36]. Later on the 
Kq(800) or k resonance in Kir scattering was also confirmed [36]. Though the mass and width of these broad 
scalar resonances predicted by different groups are rather consistent among each other, other properties about the 
resonances in the scalar family are still under a vivid debate (e.g. their nature). 

The scalar-scalar (SS) and pseudoscalar-pseudoscalar (PP) correlation functions provide us another theoretical 
tool to further study the scalar resonances [37], once the phase shifts and inelasticities in meson- meson scattering 
are well reproduced [25, 26, 27, 29, 30, 31, 32, 33, 34, 35, 38]. As derived in Refs. [39, 12, 13, 40, 20] these 
correlator functions should fulfill a set of spectral-function sum rules that imply tight constraints to the scalar and 
pseudoscalar spectra, involving nontrivial relations between them. 



We take, and further study, the meson-meson scattering amplitudes of Ref. [1] that are calculated in one-loop 
U (3) xPT plus the explicit scalar and vector resonance exchanges at tree level within the framework of resonance 
chiral theory (R^T) [41]. As a novelty, we include the contributions from the tree level exchange of pseudoscalar 
resonances in this work so that a new fit to data is also discussed. The unitarized U(3) xPT amplitudes are more 
appropriate to incorporate the heavier /o scalar resonances since the channels involving the 77 and 77' mesons can 
be taken into account (while SU(3) xPT [4] only contains the pure octet 77s held). As noticed in Ref. [26], further 
developed later in Refs. [42, 43, 44, 45, 46, 47, 48, 34, 1], the study of the resonance pole trajectories with increasing 
number of colors of QCD, Nc, is of interest to discern possible natures of the scalar resonances. E.g. a qq resonance 
evolves with a mass M <~ 0(N C ), while its width decreases as T ~ 0{Nq 1 ). For the case of a glueball its mass 
behaves also as M~ 0(N C ) with a width that decreases faster with Nc as T ~ 0(N^ 2 ) [49, 50]. In this respect 
U(3) xPT is a ls° more adequate than the SU(3) version to discuss the large Nc limit because U(3) xPT includes 
not only the pseudo-Goldstone octet, as SU(3) does, but also the singlet 771 that in the chiral limit is the ninth 
Goldstone boson at large Nc [51]. In this way, U(3) xPT has the appropriate low energy degrees of freedom in 
the chiral limit at large Nc- Additionally, as stressed in Ref. [1], the 77 meson quickly becomes much lighter with 
increasing Nc [52]. This fact has a strong impact on the dependence with Nc of the pole trajectories of the scalar 
resonances, except for the /o(500), because of its weak couplings to the channels with 7? and 77' [1]. We show below 
the Nc pole trajectories for the different resonances involved in our study from the new fit to data. 

Another interesting property of strong interactions that can be used to restrict meson- meson scattering is average 
or semilocal duality [53] . It allows one to relate the crossed channel dynamics with the s-channel one, which supplies 
another stringent constraint on the resonance properties involved in scattering. In a recent work [54] , Pelaez et al. 
have investigated how semilocal duality between the light resonances of QCD and Regge theory can be satisfied in 
7T7T scattering. These authors employ the Inverse Amplitude Method (IAM) [55, 56, 57, 58] to unitarize the SU (3) 
one-loop and SU(2) two-loop xPT amplitudes. The main conclusion of Ref. [54] is that in the physical case, i.e. 
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when Nc = 3, the / (500) resonance, with its pole position around 440 — i 250 MeV, plays a crucial role to oppose 
the vector strength from the p(770) resonance, which guarantees the fulfillment of semilocal duality. 

Having introduced the main research topics considered here for the sake of clarity we summarize our work 
and briefly comment on the main points and results obtained by considering these topics. A complete one-loop 
calculation of the scalar and pseudoscalar form factors within U (3) unitary xPT, including the tree-level exchange 
of scalar, pseudoscalar and vector resonances, is undertaken. The spectral functions of the SS and PP two-point 
correlators are constructed by using the resulting scalar and pseudoscalar form factors, respectively, which are 
unitarized for the case of the scalar ones. After updating the fit in Rcf. [1], which is also extended by including the 
explicit exchange of pseudoscalar resonances, we study the resonance spectroscopy, quadratic pion scalar radius, 
the fulfillment of semilocal duality and the spectral function sum rules in the SS — SS, PP — PP and SS — PP 
cases, which are remarkably well satisfied simultaneously. We show that it is important to take under consideration 
the high energy constraint for the coupling of the vector resonances to pseudo-Goldstone bosons, in order to keep 
semilocal duality when varying Nc- An interesting interplay between different resonances when studying the Nc 
evolution of semilocal duality and the spectral function sum rules is revealed. In the former case the scalar and 
vector spectra appear tightly related and in the latter one the same can be stated for the scalar and pseudoscalar 
ones. This issue is closely connected with the evolution of the resonance pole position with varying Nc of the 
pseudoscalar, vector and scalar resonances. 

In this respect, for larger Nc, the authors of Ref. [54] revealed that a qq subdominant component for the /o(500), 
with a mass around 1 GeV and obtained for some sets of LECs in both one-loop and two-loop calculations, is needed 
in order to fulfill semilocal duality. For those solutions, the /o(500) pole moves away from the 400-700 MeV region 
of the real axis, but it turns around moving back towards the real axis above 1 GeV as Nc becomes larger than 8 
or so. The interpretation of these results was studied preliminarily in Ref. [59], where the /o(500) is considered as 
a combination of several states and it was found that the weight of the subdominant qq component in the /o(500) 
grows with increasing Nc- Immediately, one can ask a question: apart from the /o(500) resonance, which is only 
marginally contributed by this qq component, are there any other more obvious effects from the qq seed at Nc = 3? 
If there are, what kind of role do these effects play in the fulfillment of semilocal duality? 

Concerning the first question, we provide here, confirming Refs. [26, 1], a scenario for an affirmative answer. A 
bare singlet scalar resonance Si, with its bare mass around 1 GeV, is found significant [1] to reproduce the ww 
and K K scattering data around the /o(980) energy region and, thus, becomes an important part of the / (980) 
resonance in the physical case, i.e. at Nc = 3. When increasing Nc, the physical /o(980) resonance, gradually 
evolves to the singlet Si, which acts like the role of the qq seed described in Refs. [54, 59]. Apparently, the scenario 
that we provide now is different from that in Ref. [54] , since in the latter reference the signal of the /o (980) resonance 
gradually disappears when increasing Nc- 1 Regarding the second question raised in the previous paragraph we 
find here that its role for fulfilling the requirements of semilocal duality becomes more important with increasing 
Nc, being only of little importance at Nc = 3. This behavior is in qualitative agreement with the results of [54]. 

In addition to the results presented in our short letter [37], we show our theoretical calculations in detail and 
discuss more phenomenological materials in this work, such as the fit results, resonance pole positions and their 
coupling strengths, Nc trajectories for all the relevant resonances and so on. The contents of this article are 
organized as follows. Section 2 is devoted to the introduction of the relevant chiral Lagrangian. In Sections 3 and 
4, we present the details of the calculations of the spectral functions and the quantities to measure the degree 
of fulfillment of semilocal duality, respectively. The phenomenological discussions are given in Section 5, which 
include the fit of unknown parameters to the experimental data and its consequences on the form factors, spectral 
functions, spectral-function sum rules and semilocal duality. Section 6 is then devoted to the study of the Nc 
evolution of these quantities. Finally, we conclude in Section 7. 

2 Chiral Lagrangian 

The theoretical framework that we use in the present work is U(3) xPT [60, 61], the tree level resonance exchanges 
from chiral invariant Lagrangians [41] and unitary xPT [26, 62]. The ?7(3)®f/(3) chiral symmetry in u, d, s massless 
QCD is broken because of quantum effects that violate the conservation of the singlet axial-vector current by the 
£/a(1) anomaly [63, 64, 65]. As a result the pseudoscalar rj' is not a pseudo-Goldstone boson [66]. Nevertheless, 

1 See the left panel of Fig. 1 in Ref. [54] . It is verified that the peak around 1 GeV 2 finally disappears with larger values of Nc ■ 
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from the large Nc QCD point of view the quark loop responsible for the Ua(1) anomaly [63] is 1/Nc suppressed. 
In this way, the 7/ becomes the ninth pseudo-Goldstone boson in the large Nq limit [51]. In large Nc QCD the 
singlet field 771 can be conveniently incorporated into the effective field theory by enlarging the number of degrees 
of freedom of the theory from an octet to a nonet of pseudo-Goldstone bosons, which is usually called U (3) xPT 
[60]. In this theory there are three expansion parameters: momentum, quark masses and 1/Nc, giving rise to a 
joint triple expansion 5 <~ p 2 ~ m q ~ 1/Nc- We briefly recapitulate the relevant chiral Lagrangians for our work. 
The leading order Lagrangian in U(3) xPT reads [60] 

C x = ^<u M u") + ^<X+) + ^M 2 ln 2 detu, (1) 

where (...) stands for the trace in flavor space and the last operator is the Ua{1) anomaly term that gives rise to 
the singlet t]i mass. The definitions for the chiral building blocks are 

Up = iu + DfJJu + , 
X± = u + xu + ±ux + u, 



,2 



- %/2£ 



U = u~ = e" 
D^U = d„U -i rfl U + iUlp, 

X = 2B(s + ip) , (2) 



where , , s , p are external sources [3] and the pseudo-Goldstone bosons are collected in the matrix 

' 73" 1 Te'' 8 T 7T 



$ = 



+ + ^ ^ K+ 



(3) 



The axial decay constant of the pseudo-Goldstone bosons in the simultaneous chiral and large Nc limit is denoted 
by F. In the same limit the parameter B is related to the quark condensate through (0|^g J |0) = —F 2 BS l: >. We 
do not consider isospin (I) breaking effects throughout. 

In the present work, we exploit the assumption on the resonance saturation of the LECs [41, 67], so that, 
instead of local chiral terms contributing to meson-meson scattering, we take the tree level exchanges of the scalar, 
pseudoscalar and vector resonances. In this way we keep all local contributions to meson-meson scattering up to 
and including 0(<5 3 ), while also generating higher order ones. We also include the one- loop contributions that 
count one order higher in S, as calculated in Ref. [1]. 

The terms that describe the interactions between scalar resonances and pseudo-Goldstone bosons in RxT 
are 2 [41] 

C s = CdiSsu^u 11 ) +c m (S 8 x+) +c d S 1 (u^) +c m S 1 ( X +) , (4) 
and for the vector resonances it reads 

£v = ^<VKX]>, (5) 

where the anti-symmetric tensor formalism is used to describe the vector resonances [41]. In addition, we also 
include the pseudoscalar resonances in this work, which are not considered in Ref. [1]. The relevant Lagrangian 
reads [41] 

Cp = id m (PsX-)+id m Pi(X-)- (6) 



2 The terms Cd(5gu M ){u^) + c m S\ In 2 detu were also introduced in Ref. [1] but found phenomenologically irrelevant. This is in 
agreement with the fact that from the exchange of scalar resonances they start to generate tree level meson-meson contributions that 
are of higher order, 0(S 4 ). 
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The corresponding kinetic terms for the resonance fields are [41] 

£L - -\{ v A v v v,^ - UaIv^) , (7) 

^i^^V^sV^s-M^iil) + l(^iJ 1 a M i2 1 -M| 1 i^) , (8) 

where the symbol R refers to either scalar or pseudoscalar resonances. 

One should notice that the operators appearing in Eq. (6) cause mixing terms proportional to the quark masses 
between the pseudoscalar resonances and the pseudo-Goldstone bosons. Though it is not a problem to consider 
this effect in the calculation of Feynman diagrams, it would be convenient to eliminate the mixing terms at the 
Lagrangian level. Indeed this can be accomplished by performing the following field redefinition in a chiral covariant 
way: 



dm / 1 



3x3 



Pi^Pi+i^(X-), (9) 
with J 3X 3 the unit 3x3 matrix. Substituting Eq. (9) into Eqs. (6) and (8), we have the new Lagrangian 

% = \{ V"W - M k~P~l ) + l^P^Pi - M 2 P P\) 
+ i^(V,Ps V"x- > + i^" V M Pi V"< X- ) 

-^(x-x-) + (^-^)(x-Xx-) + -, do) 

where the omitted terms, represented by the dots in the last line, denote the local chiral operators that describe 
the interactions between the pseudo-Goldstone bosons at 0(p 6 ) that we disregard throughout. We point out that 
the above procedure only eliminates the mixing terms between the pseudoscalar resonances and pseudo-Goldstone 
bosons that are linearly proportional to quark masses. The new mixing terms with higher power of quark masses 
that arise through the operators in the second line of Eq. (10) are 0(p 6 ) contributions, which are then not considered 
in this work. It is interesting to point out that after the field redefinition of Eq. (9), two new operators describing 
interactions between pseudo-Goldstone bosons appear, i.e. the ones in the third line of Eq. (10). They coincide 
with the 0(p 4 ) terms that result from the integration of the pseudoscalar resonances [41]. 

We remind that the nature of the pseudoscalar resonances is still a somewhat debated issue and also that typically 
the parameters of those resonances are not determined accurately [36] . In Ref. [68] the pseudoscalar resonances 
are generated dynamically due to the interactions between the scalar resonances and the pseudo-Goldstone bosons, 
instead of introducing them as basic degrees of freedom at the Lagrangian level. In order to compensate for 
the uncertainties arising from the not well settled properties of the pseudoscalar resonances, as well as from the 
simple tree level exchanges with bare propagators that we take here to describe them (neglecting more involved 
contributions as those pointed out in Ref. [68]), we introduce an ig-like operator [4] in our study 

^tt(x+x+ + x-x-)- (11) 

The reason behind is that only L§ an d L-j could receive contributions at the 0{p A ) level after the integration of 
the pseudoscalar resonances from Eqs. (6) and (8). Nevertheless the L7 term vanishes if one further imposes the 
large Nq relations for the pseudoscalar resonances (as we take here) 

d m = ^§ , (12) 
M Pl = M Pg . (13) 
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We stress that 5L S is different from L 8 in the xPT Lagrangian [4] and their relation can be written as 

Lf T = i^onanccs + _ ^ 

So that we interpret 5L% as the contributions from some remnant pieces that are sub- leading in 1/Nc, minding 
that the leading Nc contributions are already included in the resonance part. 

Finally, there are two relevant terms at 0(S) that only incorporate pseudo-Goldstone bosons [69] and cannot be 
generated from the exchange of the explicit resonance fields discussed above. These contributions are: 

ip — — i lndct U , 

= - 2(a^) , (15) 

with a p = (r p — Z M )/2. As we commented in Ref. [1] the Ai parameter only affects the calculation of the masses, 
scattering amplitudes and also the scalar and pseudoscalar form factors in an indirect way, i.e. through the 
renormalization of the r\\ field and its effect in the global fit is tiny. We explicitly check that if we include Ai in 
our discussion, the fitted value for this parameter approaches to zero. So that this operator is not considered in 
the following. 

For the remaining definitions of the basic chiral building blocks, the reader is referred to Ref.[l] and references 
therein for further details. 



3 Spectral Functions and Form Factors 

The two-point correlation function is defined as 

5 ab Il R (p 2 ) = i J d 4 xe lp - x < 0\T[R a {x)R b (0)}\0 >, (16) 

where the scalar and pseudoscalar densities correspond to R a = S a = q\ a q and R a = P a = iq"f^\ a q, in order, with 
A Q the Gell-Mann matrices for a = 1, . . . , 8 and A = 73x3\/2/3 for a = 0. The spectral- function sum rule in the 
chiral limit can be then represented as 

/>oo 

/ [lmU R (s)~lmU Rt (s)]ds = 0, (17) 
Jo 

where R and R' are different scalar or pseudoscalar densities mentioned above. The imaginary part of the two-point 
correlation function lmU R , also called the spectral function, is one of the key quantities that we calculate in this 
work. 



3.1 Scalar sector 

The scalar form factor of a pseudo-Goldstone boson pair PQ is defined as 

F£ Q (s) = ±(0\q\ a q\PQ). (18) 

In the present work, we focus on the chiral dynamics for the components with a = 0, 1, 2, 3, 8, which preserve the 
strangeness. By imposing the isospin symmetry, only three of the five components are independent. We take the 
neutral ones, i.e. a = 0,8,3, which correspond to the singlet and 7 = 0, 1 octet SU(3) densities, respectively. The 
scalar form factors with components a = 4, 5, 6, 7 correspond to the strangeness changing ones with I = \ and |, 
that were studied in [23] . 

The scalar spectral function is related to the scalar form factors through 

ImII s .( S ) =]>>(*) \F?{s)\ 2 6{s-st), (19) 
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where 6{x) is the Heaviside step function and the sum on i extends over the different pseudo-Goldstone boson 
pairs. In addition, s is the energy squared in the center of mass frame, s- h = (m^ + mg) 2 is the threshold of the 
ith channel and tua, tub are the masses of the corresponding two particles. In the previous equation only two- 
body intermediate states are considered, the same ones as taken in Ref. [1] to study meson-meson scattering. We 
point out that the unitarized scalar form factors, among others, include the contributions of the single resonance 
exchanges to the spectral functions. We proof below, when discussing the evolution with Nc, that the two-point 
correlator in our analysis reduces to the single resonance exchange diagram at large Nc [70]. 
The phase space factor for the i t h channel in Eq. (19) is 



„ / s _ VI s - ( m A + m B ) 2 ] [s - (m A - m B ) 2 } _ q % . . 

1D7T s o7ry s 

being qi the three momentum in center of mass frame. 

For the isoscalar case, there are five two-particle intermediate states made by a pseudo-Goldstone pair in U (3) 
%PT, namely, irir, KK, rjr), rjrf and rfrf. For the isovector case, there are three channels: irq, KK and ttt]' . The 
two-particle states with well defined isospin for the isoscalar case read 



|t"t)/=o 
\KK) I=0 



1 Itt+tt- ) + |7T~^+ ) + |7r°7T ) 

\K+K-) + \K°K°) 



V2 



\vv) 



For the isovector case they are 



\mh=o- ^ , 
|W)/=o = |W), 

W«h= = ^>. (2D 



irr]) I=1 = \irq) , 

K+K-)- \K°K°) 



\KK 



'i=i = -- 



V2 

|tiV>J=i = |W>- (22) 

We point out that the so-called unitary normalization for the identical particles, as proposed in Ref. [25], has been 
used in Eq. (21) for the 7T7t, r\r\ and r\'r( states. In this way, the normalization that we employ in this work coincides 
with the one used in [1]. We can then construct the unitarized form factors by using the partial wave scattering 
amplitudes calculated in [1] without any adjustment in the normalization. 

The perturbative calculation of the scalar form factors of the pseudo-Goldstone pairs at the one-loop level 
consists of evaluating the Feynman graphs shown in Fig. 1. As we have done for the scattering in Ref. [1], the 
contributions beyond the leading order to the form factors are calculated in terms of rj and rj' fields that result by 
diagonalizing the fields 778 and r\\ at leading order from Eq. (1) 

% = c g rj + s e r}' , 

Vi = -serj + cerj' , (23) 

with ce — cos 9 and se = sin 9. Compared to the basis with the fields r]s and r\\ , the use of the fields rj and rj' allows 
us to include the relevant Feynman diagrams conveniently in the calculation. As discussed in Ref. [1] we avoid in 
this way having to include an arbitrary number of insertions of the leading order rjs-rji mixing. As a result, one 
can consider the r\ and rj' mixing effects in the form factors order by order in the 5 counting scheme. 

There are several sources that contribute to diagram (f) in Fig. 1: rj—rj' mixing (notice that it is different from 
771-778 mixing [1], as just discussed), and the local contributions from the last two terms in Eq. (10) and the 
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(a) 





S 



< 




(d) 



(e) 



(/.) 



Figure 1: Feynman diagrams for the calculation of the scalar form factor in one- loop U(3) xPT. The wavy line denotes 
the scalar external source. The cross in diagram (d) indicates the coupling between the scalar resonance and the vacuum. 
The filled circle in diagram (f) corresponds to the vertices that only involve the pseudo-Goldstone bosons and are beyond 
leading order from Eq. (1). The different terms that contribute to diagram (f) are explained in detail in the text. 



operators in Eqs. (11) and (15). The contributions from the rj-rj' mixing to the form factors are calculated in 
the same way as for the calculation of scattering amplitudes in Ref. [1]. First, we substitute Eq. (16) of Ref. [1], 
which parameterizes the relations between 77 , 77' and the physical states 77,77', into the leading order Lagrangian 
Eq. (1). Afterwards, we calculate the form factors in terms of the physical states 77 and 77'. Due to the inclusion of 
the pseudoscalar resonances and the operator in Eq. (11), which contribute to the 77 — rf mixing parameters Si in 
Eq. (14) of Ref. [1], we present the updated results in Appendix A. In addition, the wave function renormalization 
of the pseudo-Goldstone bosons, which has been calculated in Ref. [1], also contributes to the scalar form factors. 
The explicit expressions for the scalar form factors are given in Appendix B. 

In U(3) xPT the resummation of the unitarity chiral loops is essential because of the large mass of the s 
quark and the large anomaly mass in Eq. (1). Then, in many kinematical regions the resulting pseudo-Goldstone 
thresholds are much larger than the three-momenta involved and this enhances the two-pseudo-Goldstone reducible 
loop contributions [71]. In addition, we are also interested in the chiral dynamics involving the resonance region 
where the unitarity upper bound in partial waves could be easily found to be saturated. So that it is not meaningful 
to treat unitarity in a perturbative way as in vPT. As a result, one must resum the right-hand cut that stems 
from unitarity and for that we employ Unitary vPT (UvPT), both for calculating the meson-meson partial wave 
amplitudes [26, 22, 1] and the scalar form factors [72, 73, 74]. 

The unitarization method that we use to resum the unitarity cut for the scalar form factors was developed in 
Ref. [73, 72]. It is based on the N/D method and was first applied to meson-meson scattering in Ref. [26]. For 
completeness, we recapitulate the essentials of this method here. 

In the case of two-particle intermediate states, the absorptive part of the form factor obeys the relation 



where T IJ (s) denotes the partial wave scattering amplitudes with definite isospin number I and angular momentum 
J and Z is the number of channels with the same quantum numbers IJ. The T-matrix T 00 (s), relevant to the 
scalar form factors with a = 0, 8, is a 5 x 5 matrix and T 10 (s), relevant to a — 3, is a 3 x 3 matrix. Both of them 
are studied in detail in Ref. [1] . Nevertheless due to the fact that now we also include the pseudoscalar resonances 



z 




(24) 



k=l 
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and the SL 8 operator, we need to consider their contributions in the scattering amplitudes as well. The effect from 
the pseudoscalar resonances can be easily included, since after the field redefinition the only relevant terms are 
the last two operators in Eq. (10). The pertinent expressions for the perturbative scattering amplitudes from E/(3) 
%PT including explicit exchanges of scalar and vector resonances [1] is given in Ref. [75]. The new contributions 
from the 5L 8 operator and pseudoscalar resonances (as well as from tensor resonances, introduced below) are given 
in [76], together with the expressions for the form factors. 

Following the method elaborated in Refs. [72, 73], see Ref. [74] for a simplified discussion, the unitarized scalar 
form factor can be cast as 

F I (s)=[l + N IJ (s)g IJ (s)]- 1 R I (s), (25) 

where 

N IJ {s) = T 7J (s) (2)+Res+Loop + rplj^) g IJ {g) T " (s) (2) ; 

R'is) = f / (s) ( 2 )+ Rcs + l -p + at"( s) (2) g u {s) F /( S )(2) . (26 ) 

The matrix N IJ (s) contains the crossed-channel cuts from the meson- meson scattering and the bare resonance 
poles. It was calculated in Ref. [1], and extended here by including the pseudoscalar resonance exchanges and the 
5L$ operator. The superscripts (2), Res and Loop in Eq. (26) denote the perturbative calculations from the tree 
level result using the leading order Lagrangian Eq. (1), resonance contributions (also including the operators in 
Eqs. (11) and (15)) and pseudo-Goldstone loops, respectively. On the other hand, R 1 (s) is a vector with Z rows 
constructed from the t/(3) xPT form factors similarly as N IJ (s) is calculated for scattering. The vector R 1 (s) 
does not contain any cut singularity and is real [72, 73]. 

The matrix g IJ (s) in Eq. (26) is a diagonal Z x Z matrix, with its ith non- vanishing matrix element given by 

16tt V 3 (s)i = a SL in) + log - x+ log - - x _ log , (27) 

/i z X + X- 

X± = 3 + ~ mB ± Y s \l- A < m A - » 0+ ) + ( s + m A - m l) 2 . 

where 05l(/x) is a subtraction constant. The matrix g IJ (s) collects the discontinuity caused by the two-particle 
intermediate states along the right hand cut and plays a key role in the N / D unitarization approach. We refer to 
Section IV of Ref. [1], and references therein, for detailed discussions on the calculation of the scattering T-matrix. 
The latter is given by a similar expression to Eq. (25) in terms of N IJ (s) and g IJ (s), 

T IJ {s) = [1 + N IJ {s) g IJ {s)] _1 7V 7J (s) . (28) 

Before ending this section, we introduce the scalar form factors in the quark flavor basis instead of in the singlet- 
octet flavor basis. There are two kinds of isoscalar scalar densities in the quark flavor basis: uu + dd and ss. The 
relations between the form factors in the two different bases read 

puu+dd _ _}_pa=8 _|_ pa=0 ^ (29) 

\/3 V 3 

pes = 1 pa =8 + 1 pa=0 /ggN 

y/3 V6 

The pion scalar radius (r 2 )J, an important low energy observable, is defined through the low energy Taylor 
expansion of the pion scalar form factor in the quark flavor basis [77, 78] 



F^ +dd {s) = F™+ M (0) 

with 



l + l(r 2 r s s + ... 



(31) 



mlF™ +dd (s) = 2Bm( 0\uu + dd\inT ) /=0 . (32) 

In the above equation, m stands for the up or down quark mass (we ignore isospin breaking) so that 2Bm = m 2 , 
with m n the leading order pion mass from Eq. (1). The relation between the physical mass square of the pion, 
m 2 ., and the leading order one, m 2 ., was given in the Appendix of Ref. [1]. The updated version that includes the 
pseudoscalar resonances and the 5L$ effect is collected in Appendix A. 
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3.2 Pseudoscalar sector 



For the pseudoscalar spectral function, the leading order contribution is due to the single pseudo-Goldstone boson 
exchange. The next non-vanishing contribution from the pure pseudo-Goldstone system requires at least three 
intermediate mesons, which belongs to a two-loop calculation that is beyond the scope of our current discussion. 
In order to take into account the chiral dynamics above 1 GeV in the pseudoscalar spectral function, we include 
the pseudoscalar resonance exchanges explicitly making use of the chiral invariant Lagrangian, Eq. (10). In light of 
the results of the work [68] , where the pseudoscalar resonances were dynamically generated through the scattering 
of the scalar resonances and the pseudo-Goldstone bosons, we could instead consider the scalar resonances and the 
pseudo-Goldstone bosons in the intermediate states. This is left for future work. 

The pseudoscalar spectral function is related to the pseudoscalar form factors calculated here in a simple way 

Imn P a(.s) = ^i: 5{s - m%)\H°;{s)\ 2 , (33) 

i 

where S (x) is the standard Dirac 5- function, mp i denotes the masses of the pseudo-Goldstone bosons and the pseu- 
doscalar resonances with the same isospin as the considered spectral function. Finally, Hf(s) is the corresponding 
pseudoscalar form factor 

H?(s) = 1< 0\tqj 5 X a q | Pi). (34) 

The calculation of the pseudoscalar form factor for the pseudo-Goldstone boson consists of evaluating the 
Feynman diagrams (a)-(d) in Fig. 2. Only the Feynman diagram (e) in Fig. 2 is relevant for the calculation 
of the pseudoscalar form factor for the pseudoscalar resonances. This diagram stems from the operators in the 
second line of Eq. (10). Moreover, we also consider the effects from the wave function renormalization of the 
pseudo-Goldstone bosons and express in the final results the pion decay constant F in the chiral and large Nc 
limits by the physical one F^. The latter was calculated in Ref. [1]. This reshuffling in the expressions was also 
done for the scattering amplitudes in Ref. [1]. 




(a) (6) (c) 



P 

V\* / \A= 



(d) (e) 

Figure 2: Pseudoscalar form factors of the pseudo-Goldstone bosons and pseudoscalar resonances. The symbols have the 
same meanings as those in Fig. 1, except that the wavy line stands now for the pseudoscalar external source. See the text 
for further details. 
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4 Remarks on semilocal duality 



In this section, we introduce the basic idea of Regge theory in order to formulate the notation used in this work. 
For a detailed account on Regge theory see e.g. Ref. [53]. Regge theory relates the dynamics in the i-channel with 
the s-channel for the scattering processes. In tttt scattering, the isospin well-defined amplitudes from both channels 
are related through [53] 

Tf\a,t) = \T< \s,t) +T«(M) + 

T«( M ) = + 1t«(m) - ^(M) , 

if - \TM{s,t) - ^T«(M) + ±T™(a,t) , (35) 

where the subscript of T denotes the t- or s-channel and the superscript refers to the proper isospin quantum 
number I in the t- or s-channel, respectively. 

The object of the Regge theory is the fixed-i scattering amplitude and the proposed quantity to quantify the 
fulfillment of semilocal (or average) duality between Regge theory and hadronic degrees of freedom (h.d.f.) in 
Ref. [54] is 

/ v~ n Im T t ( 5 egg >, t) dv = / v~ n Im T t ( j] adron >, t) dv , (36) 

where v = = 2s+t ~ 4 "V anc [ Sj ^ u are the standard Mandelstam variables. The "averaging" should take place 
over at least one resonance tower. Therefore, the integration region v 2 — v\ is typically taken as a multiple of 
1 GcV 2 . In this work, we shall focus on the energy region below 2 GeV 2 . Semilocal duality should be well satisfied 
for forward scattering, i.e. with t = 0, where the leading Regge trajectory (the one taken here into account) 
dominates. Like in Ref. [54], we also evaluate Eq. (36) at the threshold point t = i t h = 4m 2 to show the robustness 
of the approach under changes of t (still small compared with GeV 2 ). Different choices for n enable us to test 
the different energy regions that dominate the integrals in Eq. (36). For negative values of n, the dynamics in 
the high energy region definitely plays a more important role and it is beyond the scope of the current work. For 
large positive values of n, the integrations in Eq. (36) are then dominated by the very low energy physics, where 
resonances marginally contribute. In this work we test semilocal duality by taken n from to 3, as suggested 
in Ref. [54], which is an adequate choice for the intermediate energy region from threshold up to 2 GeV 2 , where 
several resonances contribute. 

To evaluate the left hand side of Eq. (36), the Regge asymptotic results continued down to threshold for 

Im^t Regge *) at nxe d t are needed. The explicit formulae and detailed discussions can be found in Ref. [54] and 
references therein. We do not repeat here the formulation. For the right hand side of Eq. (36), one can decompose 
the isospin amplitudes into a sum of partial waves 

ImT^M) = J2( 2J+ in™T IJ (s)Pj(z s ) , (37) 
J 

where z s = 1 + 2t/(s — 4m 2 ) is the cosine of the scattering angle in the s-channel center of mass frame and Pj(z s ) 
denotes the Legendre polynomials. By substituting Eq. (37) into Eq. (35) ImT t ; H a drons(^> ^) ^ s obtained and then 
the right hand side of Eq. (36) can be compared with the results from Regge theory. 

In order to quantify the fulfillment of semilocal duality we define, as in Ref. [54] , two types of ratios of integrals 
like these in Eq. (36), instead of comparing directly the value of the integration from the Regge theory and the 
h.d.f. . The first one, BJ n is defined as: 

r pi/- n Imr t (J) (i/,t)di/ 

K = — m • ( 38 ) 

J^u- n ImT t {I) (u,t)du 

To make closer the comparison with Ref. [54], we set in the following v\ at threshold, vi = \ GcV 2 and ^3 = 2 GeV 2 . 
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Other interesting objects to consider are the Finite Energy Sum Rule (FESR) between different isospin ampli- 
tudes with the same upper integration limit 



K 1 ' = — W— . (39) 

;; r ,-Mmif'W)^ 

=1 GcV 2 or 2 GeV 2 in the later discussions. 
Among the various cases to investigate semilocal duality between h.d.f. and Regge theory in tttt scattering, the 
golden mode is the isotensor one in the f-channel, since then the Regge exchange is highly suppressed (as there are 
no qq states with I = 2). As a result, the ratios F 21 and F 20 should tend to vanish in order to satisfy semilocal 
duality. In contrast, the dual direct s-channel allows the exchanges of several resonances with isospin I = 0, 1, see 
Eq. (35). Thus, if semilocal duality is fulfilled, one should expect the cancellation between the different resonances 
exchanged in the s-channel, which sheds light on the resonance properties. Indeed, it establishes serious relations 
between the scalar and vector spectra, as we discuss later. 

Finally, for the evaluation of the quantities in Eqs. (38) and (39) that quantify semilocal duality, the key 
ingredients are the partial waves ImT L '(s) in Eq. (37), which were studied by us in Refs. [1, 37]. This method 
is different from the one used in Ref. [54], where the resonances are regenerated by unitarizing the perturbative 
%PT amplitudes within the IAM [42, 43]. An important difference is that instead of the explicit resonance states, 
as we employ here, the IAM depends on the LECs from SU(2) or SU(3) xPT [2, 3, 4]. Nevertheless, including 
the explicit resonance states in the Lagrangian is not enough to guarantee that one can apply the theory to higher 
energies due to the important contribution from the pseudo-Goldstone boson loops. In this case, one needs to 
resum the unitarity chiral loops and a sophisticated approach based on the N/D method is formulated in Ref. [26], 
and already used to construct the unitarized scalar form factors in the previous section where we also discussed 
why U (3) xPT should be unitarized. Through the procedure to resum the chiral loops, in addition to extending the 
application energy region of RxT, one also generates physical resonances with finite widths [25, 26, 79, 80, 81, 22, 1] 
in contrast to the zero-width resonances in the bare chiral Lagrangian. In the following, we employ the partial 
wave amplitudes from this procedure to analyze the semilocal duality. 



pn« u -n ImT W^ t)dly 



5 Phenomenological Discussions 

In this section, we study the phenomenological results of the several different types of spectral-function sum rules 
presented in Eq. (17) and the ratios in Eqs. (38) and (39) that quantify the fulfillment of semilocal duality. In order 
to perform the analyses, we need to provide the values for the resonance parameters, the low energy constants in 
Eqs. (11) and (15) and the subtraction constants in Eq. (27) introduced by the unitarization procedure. 

We want to stress that all the parameters in the form factors already appear in the meson-meson scattering 
amplitudes and the pseudo-Goldstone masses. Thus, once we determine them by fitting the scattering data and 
the T) and rf masses, the form factors and the spectral functions are all predictions. 

5.1 Former Fit 

We employ here the best fit in Eq. (55) of Ref. [1], which is referred as "former fit" from now on, to calculate the 
form factors and spectral functions. In Figs. 3 and 4, we show in order the results for the pion scalar form factors 
in the quark flavor basis and the spectral functions with a = 0, 3 and 8. 
The quadratic pion scalar radius defined in Eq. (31) is found to be 

{r 2 )l =0.43 ±0.01 fm 2 . (40) 

Compared to the dispersive result from Ref. [82] 

(r 2 )l =0.61 ±0.04 fm 2 , (41) 

our result is around a 30% smaller than this well accepted value [82, 83]. The main reason for getting a smaller 
value for the pion scalar radius is that our predictions for the low energy constants L4 and L 5 are quite small by 
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using the resonance parameters from the best fit in Ref. [1] 

L 4 (m = 770 MeV) = (0.03lg;g|) x 10~ 3 , 

L 5 (n = 770 MeV) = (0.26±g;^) x 10" 3 . (42) 

These two low energy constants are important in the determination of the pion scalar radius [77]. So in order 
to improve its determination, we need to increase the values of L4 and L5 from the resonance contribution. In 
Ref. [23], it is realized that around 50% of L5 is in fact contributed by a second nonet of scalar resonances. Thus 
we decide to include a second multiplet of scalar resonances in our discussion so as to achieve a bigger pion scalar 
radius. This requires us to update the fit we did in Ref. [1], which we discuss in detail in the next section. Another 
motivation to start a new fit is that we also include the pseudoscalar resonances in the present work, which turn 
out to be important in the spectral-function sum rules that we discuss below. 

(0\uu + dd\nn)j = () (0\ss\tt-k)i =0 




400 600 800 1000 1200 1400 1600 400 600 800 1000 1200 1400 1600 

Energy (MeV) Energy (MeV) 

Figure 3: The uu + dd (left panel) and ss (right panel) scalar form factors of the pion in the quark flavor basis. The (red) 
solid lines are for the current fit and the (blue) dashed ones for the former fit. 



Current Fit Former Fit 




0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 



s(GeV 2 ) s(GeV 2 ) 

Figure 4: The scalar spectral functions for the a = (solid line), 8 (dotted line) and 3 (dashed line) correlators. The left 
panel is for the current fit and the right one for the former fit. 
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5.2 Current Fit 



Since we do not include the Air channel which turns out to be important in the energy region above 1.3 GeV in the 
isoscalar scalar case [35], we fit the scattering data in the I J = 00 channel up to 1.3 GeV as done in Ref. [1]. In this 
way, the second multiplet of scalar resonances around 2 GeV [22, 23] behaves more like a background in our study 
and we cannot fit its parameters in a precise manner. Instead, we take the values for the resonance parameters of 
the second multiplet of scalar resonances from the preferred fit of Ref. [22], given in its Eq. (6.10): 3 

c'd = V& d = c' m = V3c m ~ 40 MeV , 
M s , = M s > s = 2570 MeV . (43) 



In this reference the scattering data in the I J = ^0 and I J = |0 channels were investigated in a similar approach 
as we use here up to around 2 GeV. About the pseudoscalar resonances, we take d m — V3d m = 30 MeV and 
Mp s = Mp 1 = 1350 MeV, which lie in the commonly determined regions of these parameters [41, 20, 17]. 
In Ref. [1] we found that the scalar resonance couplings from the fit perfectly obey the large Nq relations 



Cd 



Cd 



-=vl' (44) 

so that in the current fit we also impose them. For the subtraction constants we always keep the isospin constraints 
on them [1], Finally, in the updated fit we have 18 free parameters and the fit results arc 



Cd = 


(19.8tg !]) McV 


c m — 


(41.9±|"|) MeV 




(1397±jf) MeV 


M Sl = 


(llOOljg) MeV 


M p = 


(801.2±|;|) MeV 


M K , = 


(910.0l5:S) McV 


G v = 


(62.lj£?) MeV 


10, nr] _ 
a SL — 


9 n+3.3 
z - u -4.5 


00 ,7T7r 

a SL — 


1 97+O.I2 
i - z '-0.12 


00, KK 
a SL 


u - yo -0.16 


±0,Ktt 
a SL 


-, -I o +0.12 
i - lz -0.17 


a SL 


-0.08i?;j£ 


h0,Kri' 


■ "1.25 ^ 23 


6L 8 = 


0.23+° f 9 x 10- 3 




(951+$) MeV 


A 2 = 


-0 S7+ - 19 

u -°' -0.19 


Af = 


(0.7618$) McV- 2 


c = 


i n5 +0 - 43 

l.UO_ 33 



(45) 



with x 2 /(degrees of freedom) = 784/(348 — 18) ~ 2.38. For the remaining subtraction constants, we impose the 
following relations in the fit 

00 , 777) _ 00 , 777/' _ 00, 7?' T]' _ oo . KK 

a SL — a SL — a SL — a SL ' 

20 , TTTT _ 00 . 7T7T 

a SL ~ a SL ' 

| , KTV _ \ , K 7T 

a SL — a SL ) 

IO.ttti' 10 , KK 00 . KK ,, a s 
a SL = a SL = a SL > ( 46 ) 

and all of the subtraction constants in the vector channels (which are barely sensitive to them while they are of 
natural size [62]) are set equal to a^ 71 ^". The parameters Af and c in Eq. (45) are introduced to describe the tttj 
distribution 

^ = ^Af\ T&^s) + cTl^{ S ) | 2 , (47) 
with the three momentum of the irrj system in the center of mass frame. 

3 The bare mass of the second nonet of scalar resonances is not either well fixed from Ref. [22]. E.g. in its fit (6.11) its bare mass 
is ~ 2 GeV. The range 2 — 2.5 GeV is then a realistic one for the bare mass of the second nonet of scalar resonances as follows from 
Ref. [22]. 
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The resulting plots from the fit in Eq. (45) are shown in Figs. 5, 6 and 7 by the solid curves, where we have 
also shown the best fit results from Ref. [1] by the dashed lines. The masses of rj and 77' from the new fit read 

m n = 536. 7^;;? MeV , m v , = 956t^ MeV , (48) 

which are clearly improved comparing with the ones from Ref. [1] . And the leading order mixing angle of 771 and 
r]8 introduced in Eq. (23) is 

6 = -(15.1± 2 2 ; 4 4 )° • (49) 

Compared to the former fit from Ref. [1], three more subtraction constants, a°g£ KR , a,g^' Kn and a s ^' Kv are 
set free in the current fit Eq. (45), while we reduce now the scalar resonance parameters by explicitly imposing the 
large Nc relations of Eq. (44). For the subtraction constants in the isoscalar scalar channel, o^' 11 and o°g£ KK arc 
compatible with the former fit within error bands. For the IJ = ^ channel, a s ^ ' Kn and Kn are quite similar 
and obey nonet symmetry approximately, while ag^' Kv is much more different. Nevertheless, one should notice 

that both ag^' Kr) and a s ®' Kri carry especially large error bars. About the resonance parameter a, its value from 
the new fit is larger by around 30% than the one from the former fit in Ref. [1]. This new value is closer to those 
determined in other works [26, 22, 84]. The figure c d = 19.1^2.1 MeV was reported in Ref. [26], c d = 23.8 MeV 
and Cd = 26 ± 7 MeV were given in Rcfs. [22, 84], respectively. For the parameter c m , its value from the new fit 
also increases around 30% compared to that from the former fit in Ref. [1]. However, the error bars accompanying 
c m are now considerably smaller. This is because in the new fit we impose the large Nc relation for the singlet 
couplings, Eq. (44). Concerning the bare masses, both for the resonances (scalar and vector) and the singlet 771, no 
appreciable changes in the new fit are seen, compared with the ones from the previous one [1]. The central value 
of A 2 is now around a 60% of that from the best fit in Eq. (55) of Ref. [1], though both determinations are affected 
by large uncertainties. For the vector resonance coupling Gy, the results from the current fit and the previous one 
in Ref. [1] perfectly agree with each other. Concerning 6Lg, which is introduced to compensate the uncertainties in 
the pseudoscalar resonance sector, it also carries large error bars. 

The quality of the current fit in Eq. (45) and the best fit in Ref. [1] is quite similar, as one can see from Figs. 5, 6 
and 7. Nevertheless, the prediction for L5 from resonance contributions using the new fit Eq. (45) is considerably 
increased, compared to the results in Eq. (42), due to the inclusion of the second scalar multiplet. The resulting 
values for L 4 and L 5 are now 

L 4 (fi = 770MeV) = (0.09^;^) x 10" 3 , 

L 6 (n = 770MeV) = (0.67±g$) x 10 ~ 3 • ( 50 ) 

They new value for L 5 agrees well with the recent determination L 5 = (0.58 ± 0.13) x 10~ 3 from the latest 0(p 6 ) 
SU(3) xPT fits of Ref. [85]. Concerning L4 the latest reference cannot pin down a precise value, giving the result 
L4 = (0.75 ± 0.75) x 10~ 3 . Our determination in Eq. (50) is compatible with the latter number given its large 
uncertainty. Related to the larger value for L5 in Eq. (50), the quadratic pion scalar radius from the new fit is also 
improved with the resulting value 

(r 2 )£=0.49±^fm 2 , (51) 

increasing around a 14% compared to the value in Eq. (40) from the former fit, Eq. (55) of Ref. [1]. The resulting 
scalar pion form factors and spectral functions are displayed in Figs. 3 and 4 respectively, together with the results 
from the former fit [1]. 

The resulting resonance pole positions on the complex plane from the new fit Eq. (45) are collected in Table 1. 
We refer to Ref. [1] for the discussions on how to perform the extrapolation from the physical sheet to the unphysical 
Riemann sheets. Around a resonance pole sr, corresponding to a resonance R, the partial wave amplitude Tj(s)i^j 
tends to 

Tj( S )^ _ 9R-+i9R-+j (52) 
s — s R 

By calculating the residue of the resonance pole we then obtain the product of the couplings to the corresponding 
decay modes, gn^i gn^j- The pole positions for the vector resonances p(770), if * (892) and 0(1020) agree perfectly 
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with those in Ref. [1]. Only slight changes are observed for the /o(500), /o(980), lfo(800) and ao(980) resonances. 
While all of the excited scalar resonances, such as the /o(1370), ^(1430) and ao(1450), have larger widths in the 
new fit Eq. (45) than in the previous one [1]. 
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Figure 5: Plots for the IJ — 00 case from the current and former fits. From top to bottom and left to right: the mr — > tvtt 
phase shifts (5°Sr-nnr)> the modulus of the nn — s> -kit S-matrix element (\S°% one half of the modulus of the tvtt — ¥ KK 
S-matrix element (\S^_^ K ^\/2) and the tvtt — > KK phase shifts (,S^_y K j^). The solid (red) line corresponds to the current 
fit, Eq. (45), and the dashed (blue) line represents the former fit of Ref. [1]. The error bands are represented by the shadowed 
areas, which are calculated using Eq. (45). The data for C_,„ are from Refs. [86] (green triangle), [87] (blue square) and 
the average data from Refs. [88, 89, 90] (black circle), as used in Ref. [26]. The data for |5^_,.^^| are from Ref. [88] while 
those for \S°'^_ fK ^\/2 are from Refs. [91] (blue square) and [92] (black circle). The phase shifts S°° n _^ K ^ correspond to the 
data from Refs. [91] (blue square) and [93] (black circle). 
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7l±l 
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Table 1: Pole positions for the different resonances in y/s = (M,-i§). The mass (M) and the half width (r/2) are given 
in units of MeV. The modulus of the square root of a residue is given in units of GeV, which corresponds to the coupling of 
the resonance with the first channel (specified inside the parentheses). The last two columns are the ratios of the coupling 
strengths of the same resonance to the remaining channels with respect to the first one. Note that the residues for nn, r\r\ 
and rfrf are given in the unitary normalization, see Ref. [1]. 
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5.3 Phenomenological results of semilocal duality 



We show in Table 2 the values of the ratios defined in Eqs. (38) and (39) from the current fit, Eq. (45). The 
dependence of these ratios with Nc > 3 is discussed later. We also consider the contributions of the D-waves to 
the integrals in Eqs. (38) and (39), in addition to the S- and P-waves. For including the tensor resonances, which 
originate the P-waves, we follow the formalism of Ref. [105] and take for the couplings the values determined there. 
The bare mass is adjusted such that the physical mass of the /2(1270) from the pole position agrees with the value 
in the PDG [36]. To avoid interrupting the current discussion, we give the expressions for the tensor contributions 
to meson- meson scattering in Appendix C. 

In the leftmost column in Table 2 we indicate the partial waves involved in the evaluation of the integrals in 
Eqs. (38) and (39). The values of n considered are given in the second column. In the rest of columns we give 
the quantities and P 21 , as indicated in the top row. The values t = and t = 4to 2 are used in order to show 
the stability of the results under changes in t that are small compared with GeV 2 . For the different quantities one 
should compare the numbers from Regge exchange and those obtained by including only the (S + P)-waves or in 
addition including as well the D- waves. Our results in Table 2 quantitatively confirm the conclusions of Ref. [54], 
that semilocal duality for n = 3 with 7 t = and 1 can be perfectly satisfied by including only the S- and P-waves, 
while the fulfillment for n — 2 is already marginal. For smaller values of n, the higher partial waves and higher 
cut-offs are crucial in order to satisfy semilocal duality. In this respect, we observe that once the P-waves are 
included semilocal duality is satisfied better for all the n values discussed, but particularly for n = and 1. Wc 
have also considered the role of the p(1450) but it is negligible if one takes the tttt branching decay ratio from the 
PDG [36], which is only 6%. 

For I t = 2 the situation is somewhat different. Before discussing the different numbers for P 21 in Table 2, let 
us first comment on some specific values for the ratio P 21 in order to set up a criteria that allows one to consider 
a value small and then acceptable for satisfying semilocal duality. From Eq. (35), one has that P 21 — > —1, if the 
scalar contribution is dropped (the absorptive part of the 1 = 2 channel should be negligible compared with that 
of the scalar and vector channels). In contrast, P 21 — > 1 results by neglecting the vector contribution. As we 
commented before, the ratio of P 21 should vanish if semilocal duality works well. Taking this in mind we then see 
that with the S- and P-waves, we do not find any significant signal that semilocal duality is better satisfied for a 
specific value of n, even in some cases duality is satisfied worse for a larger value of n, in contrast to the situations 
with I t = and J t = 1. However, in all cases the numbers are much smaller than 1 in absolute value, so that 
semilocal duality seems to be fulfilled quite accurately. In our scattering amplitudes higher scalar resonances are 
generated, instead of only the /o(500) as in Ref. [54], which leads to an improvement for the p — 2 channel by 
comparing the numbers for P 21 in Table 2 with the ones in Table VI of Ref. [54]. The masses of the heavier scalar 
resonances in our scattering amplitudes are close to or larger than 1 GeV 2 , as shown in Table 1. Hence, only their 
effects can be taken into account in the discussion of semilocal duality when the integration upper limit ^ max in 
Eq. (39) is larger than 1 GeV 2 . Indeed, had we set instead ^ ma x=l GeV 2 the fulfillment of semilocal duality would 
be much worse than for the ^ m ax=2 GeV 2 case, especially for n = and 1. Then, in the later discussions, we only 
consider the ratio P 21 in Eq. (39) calculated at ^ max =2 GeV 2 . On the other hand, we find that the introduction 
of the P-waves, instead of narrowing the gap between the Regge prediction and the h.d.f., worsens the situation 
for P = 2 in the n = case. It is then advisable to focus in this work on n > for P 21 [106]. 

The ratios P 20 are smaller in absolute value than P 21 because the coefficient multiplying in Eq. (35) is 
larger by a factor of two for r/ ' than for In Table 2 we do not display their values since they can not reveal 

any new information comparing with P 2 . 



5.4 Study of spectral-function sum rules 

After fixing the unknown parameters through the fit to data, we are ready to investigate the spectral-function sum 
rules presented in Eq. (17). To study them one has to include not only nonperturbative QCD dynamics but also 
perturbative QCD and operator product expansion (OPE) [5, 6]. In this way we split the integral in two parts 



/•so roo 

/ [Imn fl (s)-Imn ff (s)](is+ / [lmU R (s)-ImU R ,(s)]ds = 0. 

JO Js 



(53) 
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^max — 


2 GeV 2 


Regge 





0.225 


0.233 


0.325 


0.353 


~ 


~ 




1 


0.425 


0.452 


0.578 


0.642 


~ 


~ 




2 


0.705 


0.765 


0.839 


0.908 


~ 


~ 




3 


0.916 


0.958 


0.966 


0.990 


~ 


~ 


Ours 





0.669 


0.628 


0.836 


0.817 


-0.113 


0.040 


S + P 


1 


0.837 


0.812 


0.919 


0.908 


-0.230 


-0.087 


Waves 


2 


0.934 


0.924 


0.966 


0.962 


-0.129 


0.028 




3 


0.979 


0.976 


0.989 


0.988 


0.169 


0.345 


Ours 





0.410 


0.400 


0.453 


0.468 


0.531 


0.587 


S + P + D 


1 


0.653 


0.643 


0.694 


0.706 


0.154 


0.236 


Waves 


2 


0.850 


0.844 


0.875 


0.882 


0.027 


0.155 




3 


0.954 


0.953 


0.965 


0.968 


0.225 


0.388 



Table 2: Current Fit: R„ and are defined in Eqs. (38) and (39), respectively. In the first column from the left the 
amplitudes involved in their evaluation are shown. The different values of n are given in the second column. The rest of the 
columns correspond to R ! n and F^ 1 as indicated. Two values of t, and t t h = 4m%, are considered, as shown in the second 
row. 



The first integral, that extends along the lower-energy regime, comprises the nonperturbative region and we use our 
results in terms of h.d.f. to evaluate it. For the second one, higher in energy, the results from OPE are employed to 
evaluate the theoretical spectral functions. According to the OPE study of Rcf. [10] the different spectral functions 
with R = S, P and R' = S, P are equal in the asymptotic region in the chiral limit. 4 As a result, the second 
integral in Eq. (53) is zero. Then, testing how well a spectral-function sum rule is satisfied reduces to evaluating 
the first integral, which extends along the energy region below y/so- This is exactly the key object of our current 
study in this section. 

As discussed above, we consider the strangeness conserving scalar and pseudoscalar spectral functions for a = 
0,8,3. Hence there are 15 types of nontrivial spectral-function sum rules as those in Eq. (53). In order to 
show the results in a compact way, we display the individual values for the integration up to s for each of the 
phenomenological spectral functions in Table 3, instead of the differences between the different spectral functions. 
Note that the second integral for s > so is divergent, unless the difference between the spectral functions is taken 
as in Eq. (53). This divergent behavior is not an issue for the first integral because, as shown in Fig. 4, the 
phenomenological spectral functions are already very small for s > 2.5 GeV 2 . In this way, the results from the 
integration do not depend so much on sq as soon as they are larger than ~ 2.5 GeV 2 . This vanishing behavior 
should be expected from Eq. (19). In the latter only a finite number of two-body channels are considered so that 
if the form factors vanish for s — > oo (as expected from QCD counting rules [107]) so their contribution to the 
spectral function does. Note also that Pi{s), given in Eq. (20), tends to constant for s — > oo. The definitions of the 

4 The calculation in Ref. [10] is done up to 0(a s ) and including up to dimension 5 operators. 
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different quantities in Table 3 are 



/•so 

Wt^mir / lmHi(s)ds, (54) 
Jo 

i 

°w = E i 7 W)2 ■ ^S 8 ,S ,S 3 ,P ,P 8 ,P 3 , (56) 

where we take three different values of so to evaluate the integrations in order to show the dependences of the 
integrated results on so- The relative variance aw /W serves as a parameter to quantify how well the spectral- 
function sum rules in the scalar and pseudoscalar cases hold. 

Two situations by taking different masses for the pseudo-Goldstone bosons are investigated: the physical case and 
the chiral limit. In order to study the results at the chiral limit, it is necessary to perform the chiral extrapolation 
of our spectral functions. Though the resonance parameters in Eqs. (4)-(8) do not depend on the quark masses, the 
subtraction constants clsl in Eq. (27) introduced through the unitarization procedure vary with them. In Ref. [108], 
it is demonstrated that in the SU(3) limit (as also in chiral limit) all the subtraction constants should be equal 
for any pseudo-Goldstone pair made of the w, K and rj S mesons. Thus we need to extrapolate the subtraction 
constants from the fit, which are not necessarily equal to each other, to a common value. We hnd that at the chiral 
limit such a value indeed exists in a reasonable region (roughly from —1 to 0), where the results of the two-point 
correlators are stable and the spectral-function sum rules are better satisfied compared to the physical situation. 
In the following, we show the typical results in this region (with a common value taken for the subtraction constant 
at the chiral limit of —0.5). 

The corresponding scalar spectral functions at the chiral limit are shown in Fig. 8. 5 It is easy to conclude from 
Table 3 that the spectral-function sum rules are much better fulfilled by the new fit in Eq. (45) than by the former 
one of [1]. The most significant changes in Wi from both fits happen for the pseudoscalar cases, with i — P°, P 8 
and P 3 , which are caused by the pseudoscalar resonances that are now included in the new fit. We obtain that 
the smallest value for the violation of the spectral-function sum rules, around 10%, takes place at chiral limit by 
using the new fit result. Nevertheless, the result for the physical case from the new fit is also quite similar, with a 
violation of around 16%. 





W s o 


W s s 


W s3 


WpO 


WpS 


Wp3 


W 


cr w 


a w /W 


Physical masses 
Current Fit 
Former Fit 


8.6 9.0 9.6 
6.1 6.1 6.1 


7.4 7.5 7.7 
6.1 6.1 6.1 


7.0 7.2 7.4 
5.8 5.9 6.1 


8.9 
1.8 


11.3 
5.0 


10.1 
5.1 


9.0 
5.0 


1.5 
1.5 


0.16 
0.31 


m q = 
Current Fit 
Former Fit 


6.9 7.0 7.1 
5.2 5.3 5.5 


6.8 7.0 7.3 
5.5 5.7 6.0 


6.6 6.8 7.0 
5.2 5.3 5.5 


5.5 
0.3 


7.4 
3.0 


7.4 
3.0 


6.9 
3.8 


0.7 
2.0 


0.10 
0.53 



Table 3: Results from the integration of the spectral functions from up to so, Eq. (54). We show three results in the 
columns W s o, W s s and W S 3 by taking three different values for so: 2.5, 3, 3.5 GeV 2 . The results for the pseudoscalar cases 
are not changed for different so, since the pseudoscalar spectral functions are just some Dirac 5-functions Eq. (33). Wi with 
i = S°,S 8 ,S' i ,P°,P 8 ,P 3 , the mean value W and aw are defined in Eqs. (54)(55)(56), which are given in units of GeV 2 in 
this table. In the last column we show the relative variance aw/W. 



5 We show for later convenience the same spectral functions with Nq = 30 in the right panel of Fig. 8. 
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N c = 3 



N c = 30 
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s(GeV 2 ) 

Figure 8: The scalar spectral functions at the chiral limit for Nc = 3 (left panel) and Nc = 30 (right panel) with the new 
fit of Eq. (45). The values a —0 (red solid line), 3 (black dashed line) and 8 (blue dotted line) are considered. The inset in 
the left panel shows the same figure with reducing scale, so that the spectral functions with a = 3 and 8 are better seen. 



In the left panel of of Fig. 8 the singlet spectral function (solid line) is clearly dominated in the low energy region 
by a peak corresponding to a pole that evolves continuously with the pseudo-Goldstone boson masses from the 
/o(500) resonance pole at the physical case. This affinity of the /o(500) resonance to the singlet scalar source is in 
agreement with the study of Ref. [109], which determined that the /o(500) meson was mostly a 5(7(3) singlet. For 
the octet spectral function with a = 3 one observes neatly another low energy peak (this is better seen in the inset 
of the left panel where the scale is changed to cover more adequately the values for a = 3 and 8). This peak is due 
to the ao(980) resonance in the chiral limit, as we have checked. The higher energy peak at around 2.3 GeV 2 for 
a = 3, 8 is caused by the bare octet of scalar resonances with a common mass Ms s — 1.4 GeV, Eq. (45), shifted to 
somewhat higher energies by interference with non-resonant dynamical contributions. Finally, the strong peak in 
the singlet spectral function at around 1.4 GeV 2 comes from several sources involving the bare singlet resonance 
Si with a mass Ms 1 — 1.1 GeV, Eq. (45), and coupled channel dynamics with the ij-q' state. The lightest scalar 
resonances in the chiral limit were studied earlier in Ref. [26] and the pole positions for the Jo (500) and ao(980) 
obtained there are in agreement with ours. 

6 Results from the extrapolation of Nc 

Through the fit to experimental data in Sec. 5.2, we get the unknown couplings that appear in the chiral La- 
grangians as well as the subtraction constants. One advantage to employ the chiral Lagrangian approach in the 
phenomenological study is that once the chiral couplings are determined from one or several sets of data, we 
completely predict the other quantities that can be calculated from the same theory. For example, the resonance 
pole positions, pion scalar radius, spectral functions, spectral-function sum rules and semilocal duality that are 
discussed previously are all predictions from the fit to the scattering data. 

Another advantage by using the chiral Lagrangian approach is to study the behavior of the various quantities 
by extrapolating in the number of colors of QCD, Nc, and then to confront with the results from large Nc 
QCD [49]. This is straightforward in the chiral Lagrangian approach once the Nc behavior of the parameters 
from the chiral Lagrangians are known [41, 4]. Moreover, U(3) xPT is more appropriate to discuss the large Nc 
running compared to SU(3) or SU(2) xPT [3, 4], since the singlet T]i, explicitly included in the Z7 (3) chiral theory, 
becomes the ninth pseudo-Goldstone boson in the large Nc and chiral limits. Notice that this relevant degree of 
freedom is not treated as a dynamical active one in SU(2) or SU(3) xPT [3, 4]. It is also worth stressing that the 
rj becomes much lighter with increasing Nc [52], as explicitly shown in Ref. [1], an effect disregarded in previous 
studies [34, 42, 43, 44, 45, 46, 47, 48]. 

Only the leading order scaling with Nc of the resonance parameters, such as couplings and masses, is known 
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without ambiguities [41]. To show how robust is our knowledge on the Nc behavior for the various quantities 
studied, the sub-leading orders for the parameters in the Nc counting are necessary and could also be important 
[44, 47, 54, 1, 34]. In Rcf. [54], this uncertainty induced by the sub- leading terms in the 1/Nc expansion of the 
LECs is estimated approximately by taking different values of the renormalization scale \i. We adopt a direct way 
to estimate the sub-leading order of 1/Nc effects, that has been used in Ref. [34]. The idea is that through the fit to 
data, one can determine the bare resonance couplings and masses from the Lagrangian, which represent their values 
at Nc = 3. Once their values at large Nc are known, we perform the most general smooth extrapolation from 
Nc = 3 to the large Nc values up to and including 1/Nc suppressed corrections. Of course, the values at large Nc 
are not accessible directly by experiment, and can be ascertained only through theoretical considerations. In the last 
decades, a great progress along this line has been achieved in the analyses of short distance constraints of two- and 
three-point Green functions, form factors, r decays and irn scattering within RxT [34, 23, 67, 110, 111, 112, 113, 46]. 

The pion decay constant F n is calculated from the one- loop U(3) xPT in Ref. [1], that includes also sub- leading 
terms in the 1/Nc expansion. Throughout we always consider both the leading and sub-leading Nc scaling for 
when varying Nc as given in Ref. [1]. The values for the fit parameters in Eq. (45) are the ones taken for Nc = 3. 
On the other hand, due to the uncertainties of the values for the resonance parameters at large Nc, we consider 
several scenarios: 

• Scenario 1: We take only the leading order running with large Nc for all the resonance parameters, starting 
with their values at Nc = 3. As discussed in more detail in Ref. [1, 41] the leading running with Nc for the 
resonance parameters and meson-meson subtraction constants, asL, is given by: 

{c d {N c ), c m (N c ), G v (N c ), d m {N c )} 

= {c d (3),c m (3),G y (3),d m (3)} x 

{M Sl (N c ), M Ss (Nc),M p (Nc),M K ,(Nc),M u (N c ), M+{N C ), 

M Pl (Nc),M Ps (Nc),a SL (N c )} 
= {M Sl (3), M Sg (3), M p (3), M K , (3), M w (3), M (3), 

M Pl (3),Mp 8 (3),a SL (3)}. (57) 

For the singlet couplings Cd, c m and d m , we take the large Nc constraints in Eqs. (44) and (12). 

About the Nc running of the subtraction constant asL, we argue that it is natural to assume its constant 
behavior at large Nc, though some sub- leading Nc corrections may exist. This is based on the fact that the 
unitarized scattering amplitude, defined in Eq. (28), is in fact the sum of a series of bubble diagrams with 
the kernel N IJ (s), since Eq. (28) can be expanded as 

T IJ (s) = N IJ (s) - N IJ (s)g IJ (s)N IJ (s) + N IJ (s)[g IJ (s)N IJ (s)} 2 + .... (58) 

Within large Nc QCD it is well known that the leading Nc behavior of a meson-meson scattering amplitude, 
Nc a , corresponds to a = —1 and it can also contain other sub-leading pieces with a = —2,-3,... [49]. 
This feature for meson- meson scattering is inherited by the construction of xPT [4] . Each single diagram in 
the geometric series expansion in powers of g IJ {s) of the unitarized amplitude in Eq. (58) should decrease 
with Nc at least as 1/Nc- Focusing on the first term in Eq. (58), i.e. the kernel N IJ (s), it represents 
the perturbative results calculated from xPT and hence regardless of the resummation it should inherit the 
Nc behavior of meson- meson scattering amplitudes dictated by large Nc QCD. Its calculation within xPT 
tells us that it scales as Nq with a = — 1, including typically other sub- leading components. An immediate 
conclusion that follows is that the Nc scaling index a for g IJ (s) can be only an integer, following the above 
arguments. Moreover, a > 2 can be also simply excluded otherwise the terms with g IJ (s) in Eq. (58) could 
violate the large Nc QCD prediction to the scattering amplitudes. The case a = 0, i.e. g IJ (s) behaves as a 
constant at large Nc, is indeed the natural choice for the following reasons: 

i) The relative size between a term and the next one in the expansion of Eq. (58) is g IJ (s)N IJ (s). At 
leading order N IJ (s) behaves as p 2 /F 2 , being p 2 a typical soft external four- momentum squared attached to 
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the pseudo-Goldstone bosons. From the subtraction constant asL in Eq. (27) we then have the suppression 
factor 

as LP 2 ,. n , 
WFf ■ (59) 

For subtraction constants asL of 0(1) size, as the fitted values shown in Eq. (45), one then has the typical 
suppression for unitarity loops in xPT, given in terms of the chiral symmetry breaking scale A x pt = iirF 
[114]. In order to keep this interpretation with running Nc it is necessary that every unitarity loop is 
suppressed by an extra power of 1/Nc and, for that, the subtraction constants asL should be 0(N C ). 

ii) The combination 

asL (p) - log p 2 (60) 

in Eq. (27) is independent of the renormalization scale p. Let us consider another value p! for which cisl(p') = 
0. From the previous equation it follows that 

p! = pe~ asL ^l 2 . (61) 

From here it is obvious that if \ash{p)\ is too different from 1 then pi exponentially diverges for cisl(p) <C — 1 
or tends to for cisl 3> 1. In both cases one has too different values from the typical one for a renormalization 
scale in xPT, p <~ 0.5— 1 GcV, of the similar size to the previously introduced chiral symmetry breaking scale 
~ 1 GcV or the mass of the p resonance. The fitted values for the subtraction constants in Eq. (45) have 
the right size so as to keep an adequate value for pi which does not scale with Nc (nor should the clsl{p) so 
that Eq. (61) is meaningful). 

As commented in Eq. (14), since we explicitly include the resonance contributions to L 8l that grows like Nc 
in the 1/Nc expansion [41], we consider that SL 8 is just some remnant piece sub- leading in Nc- So we take 

SL 8 (N C ) = SL 8 (3) , (62) 

throughout the following discussion. Concerning the parameters A 2 and M , their leading Nc scaling reads 
[69] 

|A 2 (A C ),M 2 (A C )} = |A 2 (3),M 2 (3)} x A. (63) 

Scenario 2: Comparing with Scenario 1, we include the sub-leading Nc scaling for the vector resonance 
parameter Gy in Eq. (5), which describes the interaction between the vector resonances and the pseudo- 
Goldstone boson pairs, e.g. the p(770)7T7r coupling. The original type of Kawarabayashi-Suzuki-Riazuddin- 
Fayyazuddin (KSRF) relation [115] predicts Gy — F/y/2. This relation was also derived from the high energy 
constraint of the pion vector form factor at tree level [67] . An updated study of the vector form factor at the 
one-loop level [112] revealed a new version for the constraint: 

a,-*. ( M) 

This modified KSRF-like relation has also been confirmed in various contexts: partial wave irir scattering [1, 
46], radiative tau decay [113] and extra-dimension model for tttt scattering [116]. The large Nc value for the 
pion decay constant in the chiral limit can be deduced from the U (3) x?T study of Rcf. [1] with the current 

fit results in Eq. (45), leading to F ~ 80y^f MeV. 

We impose the constraint for Gy given in Eq. (64) at large Nc and the extrapolation function between 
Nc — 3 and Nc — > oo with 1/Nc suppressed corrections included is: 

~Nc 
3 

Nor/ 




Gv(Nc) = G v (N c = 3) 

- Gv(Nc = 3)-G^(N c ^™) /_3_ _ \ 
Gv(N c = 3) \N C J 



(65) 
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where G^ m (N c -> oo) = G V {N C -> oo)^^ and Gy(iVc -> oo) is given by Eq. (64). Notice that 

Gy OT (Nc — > oo) is finite in the large A^c limit. With the running of as a function of Nc from Ref. [1] we 
have the numerical value 

Gy or (^Vc^oo)~46McV. (66) 
Were the sub-leading Nc scaling for Gy not considered, i.e. if 

G^(N c ^oo) = G v (N c = 3) , (67) 

then Eq. (65) reduces to the leading behavior Gv{Nc) = Gv{Nc = 3)\J as in Eq. (57). We stress that 
the extrapolation function in Eq. (65) is unique if one considers only the next-to-leading order in the 1/Nc 
scaling for the considered parameter. Similar extrapolation functions are also used for the other resonance 
parameters when needed, as specified below. For the other parameters, we keep the same setups from Scenario 
1. 

• Scenario 3: Here, in addition to Eq. (65) of Scenario 2, we also assume that the bare masses of the p(770) 
resonance and the singlet scalar resonance Si (an important component of the /o(980) at Nq = 3) approach 
to the same value at large Nc- We can realize this scenario by increasing the bare p(770) mass by a 16% 
and decrease the bare Si mass by another 16%, so that their large Nc masses meet around 930 MeV. This 
value is indeed quite close to the preferred one for the p(770) in Ref. [54] in the large Nc limit. We take as 
extrapolation function the analogous one to Eq. (65): 



M 2 {N C ) = M 2 (N C = 3) 



M 2 (jV c = 3)-M 2 (iVc7->oo) /_3_ , 
M2(N C = 3) [n c 



(68) 



with M either the p(770) or Si bare mass. 



Scenario 4: On top of the considerations in Scenario 3 we now consider the effects of the D-waves, which in- 
clude additionally the contributions from the tensor resonances. For their resonance parameters in Eq. (C.13), 
we take the leading order scaling with Nc, 

g T (N c ) = <?t(3) x M T (N C ) = M T (3) . (69) 



We summarize the different Scenarios 1-4 in Table 4. We also considered another situation in which together 
with the characteristics of Scenario 3 we take at large Nc the mass of the octet of scalar resonances to be the 
same as that of the Si and p(770). However, we checked that this new addition produces negligible contributions 
to the ratios and . The reason is because the coupling of the octet of scalar resonances to tttt is suppressed 
numerically compared with that of the singlet scalar resonance. Due to the fact that at large Nc the qq resonances 
fall down to the bare mass position in the real axis, it will cause noticeable changes for the octet resonance pole 
trajectories. Nevertheless since the reason is obvious, we do not discuss any further this scenario. 



6.1 Semilocal duality for N c > 3 

For all the scenarios we plot in Fig. 9 the Nc trajectories of the ratio F 21 with t = 4m 2 , defined in Eq. (39). We 
verify that the results for t = are quite similar. The (red) solid line is for Scenario 1, the (green) dashed line is for 
Scenario 2, the (blue) dot-dashed line corresponds to Scenario 3 and the (magenta) dotted line does to Scenario 4. 
Between the first three scenarios the best is the third one since then the curves have the smallest absolute values 
for most of the Nc axis, as required by the Regge theory and semilocal duality. The only exception is F 21 where 
Scenario 2 gives smaller values. Concerning Scenario 4, though it gives better results for n = 2, 3 than the others, 
it leads to too large values for the n = case. In Fig. 10, we show the imaginary part of the amplitudes with well 
defined isospin in the i-channel for Scenario 3 at Nc — 3 and Nc — 30. Their integration, according to Eq. (39), 
gives F 21 . For Nc = 3 one can clearly see in the bottom-left plot of Fig. 10 a resonant bump in the 1-3 GeV 2 
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D-wave 


Scenario 1 








Scenario 2 
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Scenario 3 
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Scenario 4 
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V 



Table 4: Description of Scenarios 1-4. In the second and third columns the symbol — denotes that the sub-leading Nc 
scaling for the corresponding parameters (indicated in the first row) is not considered. In turn, yj denotes that the sub- 
leading Nc scaling is taken into account. In the last column, the symbol — means that we do not consider the contribution 
from the D-waves and \J indicates that the latter are taken into account. 



region, which is absent in Ref. [54]. This is mainly contributed by the higher scalar resonance /o(1370), and it 
plays an important role to balance the contribution of the p(770) resonance. Nonetheless, its contribution becomes 
less important with increasing Nc and for n > 2 it only has a marginal contribution for all the Nc values. When 
Nc grows, the / (500) resonance pole obtained in unitarized U(3) RxT moves deeper and deeper in the complex 
energy plane and thus barely contributes. The role played in Ref. [54] by the subdominant qq component for the 
/o(500) with a mass around 1 GeV to cancel the p(770) contribution for n = 2, 3, is played in this work by the 
/o(980) resonance, which gradually evolves to the scalar singlet Si resonance and starts behaving as a q~q state for 
Nc > 6. In both works, a qq scalar state with a mass around 1 GeV is needed in order to satisfy local duality. The 
evolution of the resonance poles with increasing Nc will be discussed in detail in the next section. 

Focusing on the solid (red) lines in the four panels of Fig. 9, resulting from Scenario 1, one can immediately 
conclude that semilocal duality, though satisfied at Nc — 3 for all the values of n, it is not well satisfied at 
large Nc for n > 1. Indeed, the situation taking place for Scenario 1 is quite similar to the one- loop IAM case 
of Ref. [54] where the /o(500) does not show a sub-leading qq component. For n = 0, both the one- and two- 
loop IAM results satisfy approximately semilocal duality at large Nc, while only the results which include the 
qq sub- leading component for the /o(500) [54], show a clear sign of duality for n = 1,2 and 3. The appearance 
of this qq subdominant component, which approaches to the real axis in the complex energy plane at large Nc, 
plays the fundamental role in the fulfillment of duality in Ref. [54]. However, one should also keep in mind that 
the iVc-behavior for the p(770) resonance in the two-loop result, with its pole position at ^/Ip~ = M p — iT p /2 = 
(950.0 — i 34.8) MeV for Nc — 12, is quantitatively different compared with the one-loop result for the same value 
of N c , y/Sp- = (710.5 - i 18.4) MeV [54]. Due to the small ratios of T p /M p for N c = 12 (in both cases it is less 
than 0.1), a reasonable approximation consists of using the narrow resonance approximate formula to estimate the 
p(770) contribution in both cases to ImT /J=11 (s). It reads [53] 

I m T^s) = 16n*^6(s-M 2 p ), (70) 



where <j(s) = y/1 — 4m 2 /s. We have adjusted the normalization of the partial wave amplitude in Eq. (70) to the 
one defined in Ref. [1], and used in this work. Now, one can analytically calculate the p(770) contribution to the 
right hand side of Eq. (36), which finally enters in the ratios of Eqs. (38) and (39). By applying Eq. (70), the 
p(770) contribution to the FESR for J t = 2 in terms of its width and mass is 

t,m;- 2 " 



f 



is- n lmTl= 2 (v,t)dv = -24Tr 2 ■ (71) 



where the integration region between v\ and ^ max always covers the p(770) mass. For simplicity we show in the 
previous equation the result at t = and similar results can be straightforwardly deduced for other values of t. 
The ratio of the p(770) contribution between the two- and one- loop cases from [54] as follows from Eq. (71) is 



rp. two — loop ( M p . two — loop \ ^i^p. one — loop) 

Tp, one-loop V M p> one-loop / ^(-^p,two-loop) 



(72) 



Taking into account the p(770) pole positions for Nc — 12 from Ref. [54], that were explicitly shown above, the 
ratio in the previous equation is 0.43 for n = 3. This implies that the p(770) contribution is reduced by more than 
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Figure 9: F% l (t = Am^). The solid (red) lines correspond to the current fit and only the leading order of the Nc scaling 
for the resonance parameters is considered, i.e. Scenario 1. The dashed (green) lines additionally include the sub-leading 
Nc scaling for Gv, i.e. Scenario 2. The dot-dashed (blue) lines correspond to take into account the sub-leading Nc scaling 
for Gv, M p and Ms 1 , i.e. Scenario 3. The doted (magenta) lines Scenario 4 so that the D-wave contribution is included as 
well. 



50% in the two-loop result compared with the one-loop case in the IAM study of Ref. [54]. When increasing Nc 
this ratio will stay put since the p(770) resonance already starts to behave as a standard qq resonance at Nc = 12, 
with its mass approaching to a constant and its width decreasing as 1/Nq [42, 43, 54]. Though this Nc behavior 
of the p(770) pole is not essential in Ref. [54] to satisfy local duality, it definitely helps to improve it. 

This reduction of the p signal with increasing Nc is explained by our present approach in a quite transparent 
way. For that one needs to take into account the sub-leading Nc scaling of the resonance parameters. It is also the 
case that when the latter are taken into account semilocal duality is also better fulfilled. Among the estimates of the 
sub-leading scaling for various resonance parameters, Gy is the most reliable one, since it can be directly derived 
by requiring a proper high energy behavior of the partial wave amplitudes, which are the key input in the study of 
semilocal duality. In addition, this constraint has also been confirmed in different processes, as already discussed 
above [112, 1, 46, 113, 116]. It turns out that the large Nc condition Eq. (64) considerably improves the fulfillment 
of semilocal duality. Thus, we provide another hint to confirm this constraint. This improvement is displayed in 
Fig. 9 by the difference between the solid (red) lines and the dashed (green) lines. The slight readjustment of the 
bare masses for the p(770) and Si at large Nc seems to improve the situations for n = and n — 3, which can 
be seen from the differences between the dashed (green) lines and the dot-dashed (blue) ones. It causes significant 
effects for n — by alternating the sign of the ratio, nevertheless the magnitudes are always less than 0.2. 

Finally, we comment on the £>-wave effects. Though the Z?-wave amplitudes can gain contributions from many 
sources, such as chiral loops, scalar, pseudoscalar and vector resonances exchanged in crossed channels, one expects 
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n = 0, t = 4m^ 



n = 2, t = iml 




Figure 10: For Scenario 3, we show ImTj 1 '^ = 4mJ) in the two panels of the top row and \mT^ 2 \t = 4mJ) in the ones 
at the bottom. The results from Nc = 3 (red solid lines) and 30 (dashed blue lines) are displayed for all the four figures. 

from phenomenological reasons that the most important ones correspond to the tensor resonances [117, 118]. From 
the change between the dot-dashed (blue) lines and the dotted (magenta) ones in Fig. 9, one can discern the role 
that the D-waves play in the FESR. As one can see, it is significant for n = 0, 1, 2 and slight for n — 3 (and it 
should be even smaller for larger values of n since then low-energy physics is enhanced, as already commented). 
For n = 1, though the tensor resonances alternate the sign of the ratios, the magnitudes are still quite small when 
varying Nc- It also clearly improves the condition for n = 2 and n — 3, though for n = it clearly deteriorates the 
fulfillment of semilocal duality. We observe that the D- waves give rise to a large contribution that overbalances the 
one of the p(770). This hints that higher vector resonances are necessary to cancel the D-wave contributions. In 
order to give a rough idea on whether the p(1450) resonance in PDG [36] can counteract the D-wave contribution, 
we include another heavier vector resonance in 7T7t scattering following the formula in Eq. (70). We refer this 
resonance as p' in the following. The only difference now is that T p i in Eq. (70) should be understood as the 
partial decay width Tp^^^. Moreover, since we focus on semilocal duality below 2 GeV 2 , we simply set the p' 
mass as 1350 MeV in order to cover its peak for the integral in Eq. (36). We find the p(1450)-type resonance gives 
a negligible contribution for all the values of n, because of its too small decay branching ratio to tttt (around a 
6%), as follows from the present information in the PDG [36]. We verify that the decay branching ratio to tttt of 
the hypothetical resonance p' needs to be around 40% in order to decrease the large contribution by the D-wave 
for the n = case down to 0.2. This is a hint in favor for the existence of a heavy p' with a significant branching 
decay ratio to 7T7t. For n = 1,2,3 no cancellations or just small ones are in fact needed in order to have a rather 
suppressed FESR with I = 2 in the t channel. Nevertheless, if included for these values of n the resulting curves 
are equally satisfactory (even better for n = 1). 
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6.2 Resonance pole trajectories with varying Nc 



Unless the opposite is stated, all the Nc pole-trajectory evolutions studied from now on correspond to the Scenario 
3 in Table 4. This is also the scenario that satisfies best semilocal duality, as seen in Sec. 6.1. 

Due to the Nc behavior of the singlet r/i mass Mg in Eq. (63), a novel feature of E7 (3) xPT, as compared with 
the SU(3) version, is that the masses of the pseudo-Goldstone bosons at large Nc, especially for 77 and 77', can be 
significantly different from their physical values at Nc — 3, as shown in Ref. [1]. On the other hand, the masses 
of the pion and kaon barely change when varying Nc- Similar results are obtained from the new fit in Eq. (45), 
that are depicted in Fig. 11. One can see there that for Nc = 30, the 77 and rj' masses decrease from their physical 
values down to around 300 MeV and 700 MeV, respectively. We point out the 77' retains a somewhat heavy mass 
at large values of Nc mainly caused by the kaon mass [1] . This is easily seen by the leading order expression of the 
77 and ?/ masses from U (3) xPT [1] . They result from the tree level calculation using the leading order Lagrangian 
Eq. (1) at the large Nc limit, i.e. M -> 0: 



2m 2 K 



(73) 



where and m v / stand for the 77 and 77' masses, respectively. 
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Figure 11: Nc running for the masses of pseudo-Goldstone bosons. 



The expressions that relate m 2 , and m 2 K with m 2 and m\ at the one-loop level are given by Eqs. (Al), (B2), 
(B4) and (B7) of Ref. [1] to which one must add the new contribution Eq. (A.l), due to 5Lg and the exchange of 
the pseudoscalar resonances. The numerical values from the new fit in Eq. (45) for the pion and kaon leading-order 
masses extracted at the one-loop level are 

m n = 136.41?;? MeV , m K = 499. 6±j^ MeV , (74) 

which leads to 

ra, = 136.41?;? MeV , m n , = 693. 3lg j MeV , (75) 

according to Eq. (73). We point out that at the chiral and large Nc limits, all the masses of the nonet of Goldstone 
bosons vanish. The Nc evolution of the leading order mixing angle in Eq. (23), given explicitly in Eq. (B7) of 
Ref.[l], is displayed in Fig. 12. At large Nc it corresponds to ideal mixing, as it should. 
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Figure 12: iVc running for the leading order r\ — rj mixing angle 9 introduced in Eq. (23). 



In the study of the evolutions of the resonance pole positions we always take the Nc running masses for the 
pseudo-Goldstone bosons as shown in Fig. 11. The resulting Nc trajectories of the various resonance poles and 
their residues were studied with detail in the Ref. [1] . Qualitatively speaking we do not find any significant changes 
in the trajectories by using the new fit in Eq. (45) and thus we confirm the conclusions obtained in this reference: 
the /o(500), ^(SOO) and ao(980) resonances go deeper in the complex energy plane when increasing the values 
of N c - For the other resonances, / (980), /o(1370), #5(1430), a (1450), p(980), #*(892) and 0(1020), approach 
to the real axis and behave like the standard qq resonances at large Nc- Notice that the /o(980) for not so large 
values of Nc does not follow the standard qq pattern (compare with the pole trajectories for the vector resonances 
in Fig. 15). E.g. its width indeed clearly increases up to Nc — 7. This is a signal of the fact that the /o(980) has 
also a strong contribution to its nature as a KK bound state [25] . 

Two variant approximations, named as vector reduced and mimic SU(3), were studied to explore the Nc 
trajectories of the resonance pole positions in Ref. [1]. In the vector reduced case, we freeze out the full propagators 
of the vector resonances in the scattering amplitudes and only keep the leading local terms generated from them, 
which are 0{p A ) in the chiral counting (or 0(<5 3 )). The purpose for introducing this approximation is to highlight 
the difference between using the tree level resonance exchanges and the local LECs in the meson-meson scattering 
amplitudes. As discussed in Ref. [1], we find that freezing the bare scalar resonances in the amplitudes does not 
lead to significant effects for the /o(500) resonance. This is why we only discuss the vector reduced case. In the 
approximation of mimic SU(3), instead of taking the Nc running masses for the pseudo-Goldstone bosons, we 
freeze them throughout and, in addition, we do not consider any 77 — rj' mixing terms. Thus, in this case, the r\ and 
77' correspond to the octet 77s and singlet 771 respectively. The masses of 7T, K and rj are fixed at the physical values 
and the 77' mass is taken from the leading order prediction of Eq. (1), which is around 1040 MeV. This mimics the 
conditions of SU (3) xPT that we consider for comparison with our U(3) xPT results. 
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Figure 13: Nc running of the pole positions for the /o(500) (or a) resonance. We give the results for Nc from 3 to 30 with 
one unit step. Both the full results and the ones from the vector reduced approximation are displayed. The left panel shows 
y/~s^ and the right one Ser- 



in Fig 13, we explicitly show the full (circles) and vector reduced (triangles) results for the /o(500) resonance 
pole (s,r) from Nc = 3 up to Nc = 30 in one unit step (as it will be also the case for the rest of pole trajectories). 
Though different strategies to estimate the sub-leading Nc scaling for the resonance parameters have been used in 
the current work, the Nc trajectories for the /o(500) are qualitatively consistent with those found in Ref. [1]. That 
is, the /o(500) resonance from the full calculation tends to fall down to the negative real axis in the s^-complcx 
plane and moves further and further away along the nearby of the s CT negative axis as Nc increases. For the vector 
reduced case this is not the case. The outcome from the latter resembles the /o(500) pole trajectory from the 
one-loop IAM [42, 44], and the full result trajectory is also one of the different /o(500) pole trajectories [54, 119]. 
It is important to remark that both results, as well as the ones of Ref. [54], are compatible for values of Nc < 10, 
and confirm again the results obtained in [42, 43], which predict a non dominant qq behavior for the /o(500). The 
outcome from the mimic SU(3) approximation, which is not explicitly shown in the figure, is quite close to the full 
results, indicating the insensitivity of the a resonance to the rj and r)' mesons, even when quite different masses 
result for these two particles as a function of Nc- As explained in Ref. [1] this relative insensitivity is due to the 
small couplings of the /o(500) to the r/rj, r/r/' and rfrf , even though these couplings are somewhat larger in the 
new fit than in the previous one. However, the largest couplings to 7T7t and KK are almost the same as in Ref. [1]. 
This pole trajectory clearly indicates that the /o(500) resonance has no significant qq or glueball components and 
it is in agreement with its dynamical generation from the isoscalar scalar 7T7t interactions. Ref. [38] calculates 
the quadratic scalar radius of the tr, (r 2 ) & a = (0.19 ± 0.02) — i (0.06 ± 0.02) fm 2 , so that it is concluded that this 
resonance is a compact one and the two pions merge inside it. As a result, a four-quark picture is more favorable 
than the molecular description. 

On the other hand, the resulting pole trajectories of the heavier scalar resonances for the full result and the 
mimic SU(3) approximation are quite different between each other, as one can see in Fig 14. This tells us that 
the scalar resonances /o(980), /o(1370), Aq(1430) and ao(1450) are sensitive to the -q and states, to which they 
couple strongly as shown in Table 1. For the vector resonances, we show the Nc trajectories in Fig 15, where the 
different treatments of the // and if mesons barely change the f>(770) resonance and only a little the A"* (892). A 
kink structure has been found for the coupling strength between the A* (892) resonance and the Kr\ channel in the 
full result around Nc = 14 [1], which is caused by the crossover of the A* (892) mass by the Kr\ threshold. This 
is mainly due to the variation of the rj mass. We find the structure from the new fit is quite similar to the one in 
Ref. [1]. For the </>(1020) and w(780) resonances, we verify that it is also well behaved as a qq resonance at large 
Nc- The situation is even simpler in this case, since it only involves the A A channel in our approach. 6 For the 

e In the present study we are missing the important 3n channel for a realistic treatment of the o;(780) resonance. We can then only 
study its mass which can be fixed to its experimental value by tuning its bare mass, as similarly done for the other vector resonances. 
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Figure 14: Nc trajectories of the scalar resonances /o(980), /o(1370), JC"o(1430) and ao(1450). We show the results for 
Nc from 3 to 30 with one unit step. The results for both the full (circles) and the mimic SU(S) (squares) calculations are 
shown. The differences between the two sets of trajectories show the sensitivity of the corresponding resonances to the rj 
and r]' mesons. Specific symbols are assigned to the ao(1450) resonance for Nc = 4, 5, 6, 7 in order to show more clearly its 
pole trajectory. 



Kq(800) and ao(980) resonances, we re-confirm as in Ref. [1] that their pole positions go deeper and deeper in the 
complex energy plane when increasing Nc, which are displayed in Fig 16. 

6.3 Spectral-function sum rules for Nc > 3 

As a further application, we study the Nc evolution of the spectral functions and the spectral-function sum rules 
in this section. The scaling with Nc for the two-point correlators, i.e. Wi in Eq. (54), is proportional to Nc in the 
large Nc QCD [49]. The situation for the pseudoscalar correlators is obvious, since after substituting Eqs. (33), 
(B.ll) and (B.12) into Eq. (54) one has 

Wi cxF^ + 8d 2 m + --- , z = P°, P 8 , P 3 , (76) 

where the ellipsis stand both for sub-leading Nc or suppressed pieces in the chiral counting. From Eq. (57) and the 
fact that Fn scales as \/Nc, it results that the leading Nc scaling of the pseudoscalar correlators is proportional 
to N c - 

However, the situation for the scalar correlators is not so obvious. In fact, if we simply use the tree level results 
in Eqs. (B.6), (B.7) and (B.8) from the Lagrangian in Eq. (1) and substitute them into Eq. (19) to calculate the 
correlators in Eq. (54), one would conclude that the latter behave as a constant at large Nc, instead of running like 
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Figure 15: Nc running of the pole positions for the /o(770) and K* (892) resonances. For notation see Fig. 14. 
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Figure 16: Nc dependences of the Kq(800) and ao(980) resonances both in the fourth Riemann sheets from Nc = 3 to 30 
with one unit step. See Ref. [1] for details on the definitions of the different Riemann sheets. 



Nc- The subtlety, as shown in detail below, comes from the fact that at large Nc the width of a qq-like resonance 
tends to zero and the imaginary part of its propagator behaves then as a Dirac 5-function, which is the dominant 
contribution to the correlators in Eq. (54). 

In the following discussion, we demonstrate that by properly taking into account the resonance contributions 
in the unitarized scalar form factors Eq. (25) for the pseudo-Goldstone boson pairs, the scalar correlators in our 
formalism also scale like Nc, as required by large Nc QCD. In order to make the analytical discussion neatly, we 
illustrate the proof in the chiral limit and consider the single channel case, i.e. we only include the 7T7r channel 
both for the scalar form factor and the scattering. The relevant parts -R 7 (s), N IJ (s) and g IJ (s) in Eq. (25) are 
just functions not matrices for the single channel. Notice also that in the large Nc limit, only the tree level 
contributions, i.e. the (2)+Res parts in Eq. (26), can survive. We focus in the case a = 0, while the others can be 
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R a=0 (s) = 


-2- 


N 00 (s) = 


s 

F 2 " 




2c 2 M; 


+ 


F 4 s 



(77) 



obtained analogously. The corresponding expressions in the chiral and large Nc limits are: 

8c rf c m s 
F 2 Mj - s ' 

^ + k l0g[1 + uf 

" [1 + + J~ 4 J^T— S , (78) 

where we assume exact large Nc nonet symmetry for the scalar resonances in Eqs. (4) and (8), i.e. imposing that 
the masses for singlet and octet scalar resonances are equal and taking the large Nc relations in Eq. (44). In 
the previous equation Ms denotes the bare mass of the scalar nonet. In fact for the scalar resonances, only the 

combination ^J^S$ + \J~^So is the one relevant in the large Nc limit for the irir form factor and scattering. The 

interaction between the other orthogonal combination y^<So — \f^^s an d tttt is 1/Nq suppressed. 

Let us concentrate in the s region around the bare resonance pole at s = M| in Eqs. (77) and (78). Substituting 
Eqs. (77) and (78) into Eq. (25), we have the simple expression around the resonance pole s — > Mj for the unitarized 
form factor 

o _ 8c d c m M 2 1 

" 1 j ~ F 2 M 2_„_,3c^£i' 1 ' 

1V1 S S 1 16wF 4 

where we have dropped the tiny contribution that stems from the real part of the function g (s) (that would 
give rise to a self-energy contribution to the resonance bare mass), and explicitly show the contribution from the 
imaginary part in the chiral limit 

W» = -^. (80) 

By substituting the tree level decay width of S = ^J^S$ + ^/fSo to tttt calculated from the Lagrangian in Eq. (4) 
at the chiral limit [46] 

rs= 16^' (81) 



into Eq. (79), we get the standard Breit-Wigner propagator for the scalar resonance 

F ™ {S) F 2 Ml- S - i M s T s ■ {8 ) 

Due to the fact that T$ behaves as l/Nc in the large Nc limit, we can use the standard narrow width approximation 
to write the form factor squared in terms of the Dirac <5-function 

KT o (s)|2 = 1 



F4 (M| - a)2 + (M s T s y 



limK2 clS(s-Ml), (83) 



3 

where we have used Eq. (81) and the following way to approach the Dirac S- function: 

1 M S T S 



7T (MJ - s) 2 + (M s r s y 



5(s - M 2 S ) , when M S T S -> . (84) 



By combining Eqs. (19), (83), (54) and (57) we obtain the leading Nc behavior for the singlet two-point scalar 
correlator 

W s o ex c 2 „ ex N c , (85) 
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which has the same scaling as the pseudoscalar case Eq. (76). This is important because otherwise we would run 
into the contradiction that it would not make sense to consider the spectral-function sum-rules as a function of 
Nc- The generalization to the cases with a — 3 and 8 is straightforward. 

We stress that we can make a close analytical discussion only at the leading order in Nc and in the single channel 
case. Beyond these two conditions, it is rather difficult to perform the discussions analytically even in the chiral 
limit. Instead, our full calculations for the two-point correlators which include both the sub-leading Nc parts, 
such as the chiral loops, and the multiple-channel dynamics, are carried out numerically. We show Wi x Z/Nc in 
Fig. 17 for i = S°, S s , S 3 (left panel) and P°, P 8 , P 3 (right panel). The flat behavior at large N c clearly shows 
that all the two-point correlators Wi are proportional to Nc for large values of Nc- In addition, the asymptotic 
value is reached rather quickly with Nc, establishing already for Nc ^ 10. The Wi results are nearly the same at 
large values of Nc either for the SS correlators or the PP ones as it is required by the spectral-function sum rules, 
Eqs. (17) (the discussion after Eq. (53) should be taken into account). To be more precise, we find that the relative 
variance aw /W for the six points at Nc = 30 in Fig. 17 is only 5%. This indicates that the spectral-function sum 
rules arc better satisfied in large Nc than for the Nc — 3 case, reported in Table 3 with a relative variance of 10%. 

The reason behind this improvement can be attributed to the fact that at large Nc the correlators reduce to the 
single resonance or pseudo-Goldstone exchanges and the large Nc constraints from the SS — PP sector requires 
that (as previously worked out in Ref. [20]) 

<4 = <L + ^" ■ (86) 

This is perfectly satisfied by the results from our research, with c m determined from the fit in Eq. (45), d m — 30 MeV 
adopted in this work [41, 20, 17] and F ~ 80 MeV predicted from the scalar resonance parameters in the fit Eq. (45) 
using the U(3) x?T one-loop calculation of Ref. [1]. 

It is also worth considering the evolution with Nc of the scalar spectral functions. In the right panel of Fig. 8 
the latter are shown for Nc = 30. It is clear that the spectral functions are then completely dominated by the 
bare scalar pole S% and the octet of scalar resonances at around y/s = 1.4 GeV. The low energy peaks in the left 
panel of Fig. 8 disappear for large Nc due to its meson-meson dynamically generated nature, in agreement with 
what was already shown for the pole trajectories with Nc in Sec. 6.2. 

All the points in Fig. 17 are calculated by setting the upper limit of the integral in Eq. (54) to sq — 3 GeV 2 . 
We verify that the results are quite stable for s = 2.5 GeV 2 and 3.5 GeV 2 , as it was already discussed for Nc = 3 
in Table 3. 
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Figure 17: Wi x 3/Nc as a function of Nc within Scenario 3. We distinguish between the chiral limit and the physical 
case as indicated in the panels. We also show the corresponding value of a for each point. 
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7 Conclusions 



In this work, we perform a complete one-loop calculation of the strangeness conserving scalar and pseudoscalar 
form factors within U (3) xPT. We employ the resonance chiral Lagrangian to include the tree level resonance 
exchanges, instead of using the LECs from the local operators at higher chiral orders (which are supposed to be 
saturated by the resonance exchanges within the assumption of resonance saturation [41]). We extend the previous 
one-loop calculation on meson-meson scattering [1] by including the pseudoscalar resonance exchanges. The full 
results for the scalar form factors are calculated by using the unitarization method in Rcfs. [73, 72], based on 
the N/D method, taking as input the perturbative calculations of the scalar form factors and the partial wave 
scattering amplitudes. Then, we employ different form factors to calculate the various spectral functions, which 
allows us to study the spectral-function sum rules in the scalar and pseudoscalar sectors. 

The unknown parameters in our theory are fitted to the experimental data. Comparing with the discussion 
in Ref. [1], a new fit is carried out in the present work, due to the inclusion of new ingredients, such as the 
pseudoscalar resonances and the second multiplet of scalar resonances. Various poles in the complex energy plane 
for the / (500), /o(980), /o(1370), a (980), a (1450), if *(800), if *(1430), p(770), K*(892) and 0(1020) resonances 
are found and they agree well with the values from PDG [36] . The coupling strengths of the various resonances 
to the pseudo-Goldstone boson pairs are also calculated. The resonance content of the new fit presented here and 
that of Ref. [1] is quite similar, both at the qualitative and quantitative levels. However, the new fit gives a better 
numerical value for the pion scalar radius than the one in Ref. [1]. The former also contains the pseudoscalar 
resonance exchanges which implies that the spectral-function sum rules in the scalar and pseudoscalar cases are 
satisfactorily fulfilled at the level of around a 10%. In addition, the new fit results are employed to investigate 
scmilocal duality in irir scattering between Regge theory and the dynamics in terms of hadronic degree of freedom. 
We conclude that semilocal duality is well satisfied in general terms for n > 1. 

An important advantage of working within the chiral Lagrangian approach is that it allows us to study the Nc 
evolution of the various quantities calculated, once the Nc scaling behavior of the parameters in the theory is 
settled. The leading evolution with Nq is known unambiguously, while this is not the case for the sub-leading one. 
We propose three scenarios to take into account the uncertainties of the sub-leading Nc scaling for the resonance 
parameters and another one to include the tensor resonances in the study of semilocal duality. Interestingly, we find 
that semilocal duality at large Nc imposes strong constraints on the resonance parameters and Scenario 3 turns 
out to be the best one. Under this scenario we then study the extrapolation of Nc > 3 for many other quantities: 
the masses of the pseudo-Goldstone bosons, the leading order 77 — 77' mixing angle, the various resonance poles 
and the two-point correlators. We find that semilocal duality and the spectral-function sum rules in the scalar 
and pseudoscalar cases are well satisfied by using the parameters from the fit to experimental data, both for the 
physical case at Nc = 3, and for varying Nc > 3. It is important to stress that the fulfillment of semilocal duality 
and the spectral-function sum rules in the scalar and pseudoscalar sectors for both the physical and the large 
Nc cases is not trivial and gives a deep insight in the isoscalar scalar spectrum for ^/s < 1 GeV. In the physical 
case, the /o(500) resonance plays an important role for both semilocal duality and spectral- function sum rules, 
dominating the strength in the low-energy scalar isoscalar channel. However, when Nc increases this resonance 
moves deeper in the complex energy plane and its contribution tends to vanish. Instead, the /o(980) resonance, 
which gradually evolves to the singlet scalar bare resonance Si when increasing Nc, plays the key role for the 
fulfillments of semilocal duality and spectral-function sum rules. It is then clear that the physical picture for the 
scalar isoscalar sector in the real and large Nc worlds are very different. The former case is dominated by the 
/o(500) resonance, which should not have a significant qq or glueball components as it is clear from its Nc pole 
trajectory. Within our approach it is dynamically generated by the 7T7t self-interactions. On the other hand, in the 
large Nc limit the lightest resonance is the singlet Si bare or pre-existing resonance with a mass around 0.9-1.0 
GeV, which for Nc = 3 is a component of the /o(980). 
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In this section, we provide the contributions to the self-energies of the pseudo-Goldstone bosons that stem from 
the exchange of the pseudoscalar resonances, i.e. the operators in the last line of Eq. (10), and the SL 8 operator 
in Eq. (11). About the second multiplet of scalar resonances, its contributions share the same form as the lowest 
multiplet given in [1] with the obvious changes in couplings and masses. Hence, we do not reproduce them here. 
Notice that one needs to combine the results in the present section and those in the Appendices of Rcf. [1] to get 
the full expressions. 

We point out that the exchange of pseudoscalar resonances and the SLg operator do not contribute to the wave 
function renormalization of the pseudo-Goldstone bosons and the kinetic terms of r}-rf mixing parameterized by 
S v , S v ' and 5k in Eq. (14) of Rcf. [1]. For the pion and kaon masses these new contributions read 



Appendix 



A Masses and mixing parameters 




(A.l) 



(A.2) 



The mass terms in the 77-77' mixing are: 




+ 



+ Sg(4m% - 4m 2 K ml + 3m 4 ) 



(A.3) 




8d 2 m [8c 2 e (ml - m 2 K ) 2 - 4V2 C gs (Am 4 K - hm 2 K ml + m 4 ) + s 2 e (ml - Am 2 K ) 2 } 

9FZM 2 P& 

M 2 m [c 2 {ml + 2m 2 K ) 2 - A^c s e {2m 4 K - m 2 K ml - mj) + 8s 2 (ml - m 2 K f\ 

mM 2 Pl 



+ 




cl(4m 4 K - Am 2 K ml + 3m^) - 8\^cgsem 2 K (m 2 K - ml) 



+ s 2 g (8m 4 K - 8m 2 K ml + 3m 4 T ) , 



(A.4) 



38 



8d 2 



9F*M% 



2V2c 2 {4m 4 K - 5m 2 K ml + m%) + c e sg(-8m 4 K - 8m\m\ + 7m 4 ) 
2V2s 2 e {4m 4 K ~ 5m 2 K ml + ml) 



+ 



8dl 



3FZM 2 Pi 



+ 



2V2c 2 e (2m 4 K - m 2 K ml - m%) + c g s (-4m 4 K + 2§m 2 K m\ - lm%) 

+ 2V2s 2 g (-2m 4 K + m 2 K ml + m 4 .) 
64SL 8 m 2 K (ml - m 2 K )(V2c 2 - c e sg - y/2s 2 g ) 

m 



(A.5) 



where 5 m 2_, S m 2 and S m 2 are defined in Eq. (14) of Rcf. [1]. Notice that the two operators in the last line of the 

pseudoscalar resonance Lagrangian in Eq. (10) and the 5L 8 operator of Eq. (11) do not contribute to the pion 
decay constant F v . 



B Analytic expressions for scalar and pseudoscalar form factors 



In order to calculate the scalar form factors in the isospin basis of Eqs. (21) and (22), we need to evaluate for the 
isoscalar case eight form factors in the charged bases, namely: F°^ , F°f w _, F°'3 n+ , F^f +K _, F°' 8 _„, F^ 8 , F^ 8 , 

and F°; 8 ,. For the isovector case, one needs: F^, F^ +K _, F^ — and F% 0j)> . 

The leading order contributions to the scalar form factors defined in Eq. (18) stem from the Lagrangian in 
Eq. (1). The expressions for a = read: 



For a — 8 one has: 
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(B.7) 
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And the results for a = 3 are: 

3 _ 2( C g - V2 S g) 

F K+K- = ~ F k»k" = 1 ' 
3 _ 2(V2CQ + Sg) 

Ko v > ~ 7f • (B-8) 

The scalar form factors in the quark flavor basis defined in Eqs. (29) at leading order arc 

jpuu+dd puu+dd puu+dd ^ 

7r°-7T 7T+7T~ 7T~ 7T+ ' 

T^uu+dd j^uu+dd -i 

r K+K- ~ K°~K° ~ ' 

F^ +3d =^(V2c 2 -c e s e -V2s 2 e ), 

F™+ dd = \(s 2 e + 2V2c eSe + 2) , (B.9) 



and 



T7SS TTiSS TPSS f» 

T7ISS T7ISS 1 

^ K+K~ - t K a- K <) - 1 > 

F^ = l(c 2 e + 2V2c eSg + l), 

F^ = ~(V2c 2 e -Cg S g-V2s 2 e ), 

F^ ri , = ^(s 2 -2V2cese + l). (B.10) 
At the same order, the pseudoscalar form factors for the pseudo-Goldstone bosons defined in Eq. (34) are 

tf° = , = -2F W s 9 , = 2F n c e , 

H* = , ffS = 2F V eg , JT» - 2F W s 9 , 
H* =2F„, H* = 0, H? y =0, (B.ff) 

and for the pseudoscalar resonances they read 

-^Pseudoscalar Resonances 4\/6^m ; 

^Pseudoscalar Resonances = — 4\/2rf m , 

-^Pseudoscalar Resonances = — 4\/2rf m . (B.12) 

We provide the remaining expressions, such as those from the chiral loops, scalar resonances, pseudoscalar 
resonances, A 2 and 5L$, in the Mathematica code [76]. 

C Tensor resonances in meson-meson scattering 

We follow the framework proposed in Ref. [105] to include the tensor resonances in meson-meson scattering. The 
relevant Lagrangian reads 

Ct = -\{T^D^'^T^)+ aT {T^,u v })+^u v u v )+ 1 {T^ X +), (C13) 
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where the nonet of the tensor resonances are collected in the matrix 



T 



/ Q 2 _|_ fj I $1 n -\- 

V2 + V6 + V3 2 " 2 



K* 



Kf 



is*- 1^*0 2 /I I Ji. , 



(C.14) 



See Ref. [41] for the definition of the remaining chiral building blocks. Like in the vector resonance case [1], ideal 
mixing is also assumed for the tensor resonances: 



/ 2 (1270) = y|/ 2 1 + \/i/ ; 



3" 2 



/2 = y|/2 ■ (c.i5) 

In the present work, the purpose to introduce the tensor resonances is to discuss semilocal duality in tttt scattering. 
So we will set equal all the tensor masses, which is determined by the most relevant resonance /2(1270). The first 
operator in Eq. (C.13) corresponds to the kinetic term while the remaining terms describe the interactions between 
the tensor resonances and the pseudo-Goldstone boson pairs. The terms proportional to /3 and 7 do not contribute 
to the on-shell decay of the tensor resonances because the tensor field is traceless in the space-time indices [105]. 
Nevertheless, they do contribute to meson-meson scattering. It is argued in Ref. [105] that though the final result 
is independent of the choice of j3, it is convenient to set f3 — —qt to avoid the inclusion of the 0(jP) LECs in order 
to fulfill the high energy constraints for the forward tttt scattering. The 7 operator is always accompanied by quark 
masses and its effects should be much less important than the f3 term. So we will omit the former term throughout 
as done in Ref. [105]. 

Last but not least, one should guarantee the right high energy behavior of meson-meson scattering in the presence 
of tensor resonances resulting from the Lagrangian in Eq. (C.13). The high energy constraint imposed in Ref. [105] 
concerns the fulfillment of a once subtracted forward dispersion relation for 7r + 7r° elastic scattering, which turns 
out to play a crucial role to get the correct prediction for the LECs. We follow the same approach here to satisfy 
this high energy constraint. Moreover we calculate all the relevant meson-meson coupled channels for scattering 
in C/(3) xPT, not just 7T7t, since the other processes can enter through the unitarization procedure [1] that couples 
the different states with the same quantum numbers. 

The final results on meson-meson scattering contributed by the tensor resonances and the high energy constraints 
are given in the Mathematica code [76]. Due to the consideration of the scattering processes involving r\ and 77', 
additional LECs are needed to guarantee the proper short distance constraint imposed by the forward scattering. 
The short distance constraints that we find here are 



,>,u _ ~3t oSD _ 9t r sd _ _^9t_ oSd _ 9t , r u - 

Pl3 "M 2 ' Pu ~ M 2 ' Plb ~ M 2 ' Pl6 ~ 2M 2 ' (Uibj 



in addition to the ones determined in Ref. [105] 



dSD _ 9t oSD _ 9t oSD _ 2 ff' 2 

1 ' 2M 2 ' 12 " Ml ' Pz ~ Mr 



/*r = -f&> 03=i%- ( c - 17 ) 



The convention to label the LECs in Eq. (C.16) is the same as in Ref. [120], where one can also find the corresponding 
monomials multiplied by the coefficients with i = 1, 2, 3, 13, 14, 15, and 16. We checked that only after the short 
distance constraints are fulfilled, the introduction of the tensor resonances does not spoil the fit that we obtained 
in Ref. [1] and the new one in Eq. (45). 

For the value of gx, we take the result gr = 28 McV from Ref. [105]. As in the vector channels [1], due to the 
dominant role played by the tree level tensor resonances, one does not expect the subtraction constants arising 
from the unitarization procedure to play an important role in IJ = 02. We simply fix their values to the IJ = 00 
channel and checked that the results are stable under 0(1) changes in these numbers, so that the subtraction 
constants keep natural values in Eq. (45). The bare mass of the tensor resonances is adjusted to Mt — 1300 MeV, 
leading to the pole position y/s — (1275.2 — i 75.8) MeV, which is close to the values of the /2(1270) resonance 
given in the PDG [36]. 
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